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We review the collectivity and the suppression pattern of charmed mesons - produced 
in proton-nucleus and nucleus-nucleus collisions at SPS (~ 158 A-GeV) and RHIC en- 
ergies (~ 21 A-TeV) - in comparison to dynamical and thermal models. In particular, 
we examine the charmonium 'melting' and the 'comover dissociation' scenarios - im- 
plemented in a microscopic transport approach - in comparison to the available data 
from the SPS and RHIC. The analysis shows that the dynamics of c, c quarks at RHIC 
are dominated by partonic or 'pre-hadronic' interactions in the strongly coupled plasma 
stage and can neither be modeled by 'hadronic' interactions nor described appropriately 
by color screening alone. Both the 'charmonium melting' and the hadronic 'comover 
absorption and recreation model' are found, however, to be compatible with the experi- 
mental observation at SPS energies; the experimental ratio of ^/'/J/^ versus centrality 
clearly favors the 'hadronic comover' scenario. We find that the collective flow of charm 
in the purely hadronic Hadron-String Dynamics (HSD) transport appears compatible 
with the data at SPS energies, but substantially underestimates the data at top RHIC 
energies. Thus, the large elliptic flow V2 of D-mesons and the low Raa{pt) of J/\P 
seen experimentally have to be attributed to early interactions of non-hadronic degrees 
of freedom. Simultaneously, we observe that non-hadronic interactions are mandatory 
in order to describe the narrowing of the J/ty rapidity distribution from pp to central 
Au + Au collisions at the top RHIC energy of y/s = 200 GeV. We demonstrate addition- 
ally that the strong quenching of high-p^ J/Nl/'s in central Au + Au collisions indicates 
that a large fraction of final J/ty mesons is created by a coalescence mechanism close 
to the phase boundary. Throughout this review we, furthermore, provide predictions for 
charm observables from Au+Au collisions at FAIR energies of 25 - 35 A-GeV. 
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1. Introduction 

The formation of a quark-gluon plasma and its transition to interacting hadronic 
matter has motivated a large community for almost three decades According to 
current understanding, the universe in the 'Big Bang' scenario has evolved from a 
quark-gluon plasma (QGP) to color neutral hadronic states within the first second 
of its lifetime. In this context, the phase transition from partonic degrees of freedom 
(quarks and gluons) to interacting hadrons is a central topic of modern high-energy 
physics. In order to understand the dynamics and relevant scales of this transition 
laboratory experiments under controlled conditions are performed with relativistic 
nucleus-nucleus collisions. The study of nuclear matter under extremely high baryon 
density and temperature - where according to lattice quantum chromodynamics 
(QCD) - the hadronic matter transforms to a strongly interacting quark-gluon 
plasma (sQGP) - is the aim of a variety of experiments at current and future 
facilities: NA38, NA50 and NA60 at the Super-Proton-Synchrotron (SPS) 1214151 
PHENIX^ STAR, PHOBOS and BRAHMS at the Relativistic-Hcavy-Ion-Collider 
(RHIC) El ALICE at the Large-Hadron-Collider (LHC) EEl CBM and PANDA at 
the Facility for Antiproton and Ion Research (FAIR)^ NA61 at the SPS Heavy-Ion 
and Neutrino Experiment (SHINE) EE MPD at the Nuclotron-based Ion Collider 
Facility (NICA) M 

Relativistic nucleus-nucleus collisions have been studied so far at beam energies 
from 0.1 to 2 A-GeV at the SIS (Schwerlonen-Synchrotron), from 2 to 11.6 A-GeV 
at t he AG S (Alternating Gradient Synchrotron) and from 20 to 160 A-GeV at the 
gpg [12113] While part of these programs are closed now, the heavy-ion research has 
been extended at RHIC with Au+Au collisions at invariant energies yfs from ~ 20 
to 200 GeV (equivalent energies in a fixed target experiment: 0.2 to 21.3 A-TeV). 
In the near future, further insight into the physics of matter at even more extreme 
conditions will be gained at the LHC, which will reach center-of-mass energies of 
the TcV scale. Apart from LHC, the SPS successor SHINE will operate at CERN 
in order to scan the 10A-158A-GeV energy range with light and intermediate mass 
nuclei At FAIR, which is expected to start operation in 2015, collisions of gold 
nuclei from 5 A-GeV up to 35 A-GeV will be studied, thus exploring the high baryon 
density region of the nuclear matter phase diagram. At NICA it is planned to start 
the experimental program of colliding Au and/or U ions as well as polarized ligh t 
nuclei at energies up to of 5 A-GeV in 2013 (an upgrade to 9 A-GeV is foreseen 

In fact, estimates based on the Bjorken formula for the energy density 
achieved in central Au+Au collisions suggest that the critical energy density for 
the formation of a QGP is by far exceeded during a few fm/c in the initial phase 
of the collision at RHIC energies (y/s up to 200 GeV), but sufficient energy den- 
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sities (~ 0.7-1 GeV/fm 3 cf. Ref. H^j) might already be achieved at AGS energies 
of - 10 A-GeVH^llSDland tmis also for the energy regime of FAIR. The crucial 
question is, however, whether the partonic system really reaches thermal and chem- 
ical equilibrium in relativistic nucleus- nucleus collisions. Nonequilibrium models are 
needed to trace the entire collision history, so that one can study the nature of the 
transition and extract the characteristics of the partonic phase from data. 

Currently many properties of the new sQGP phase are still under debate and 
practically no dynamical concepts are available to describe the freeze-out of par- 
tons to color neutral hadrons that are subject to experimental detection. Early 
concepts of the QGP were guided by the idea of a weakly interacting system of 
partons (quarks, antiquarks and gluons), because the entropy density and energy 
density were found in lattice QCD to be close (within 20% accuracy) to the Stefan 
Boltzmann limit for a relativistic noninteracting system However, this notion 
had to be given up in the last years, since experimental observations at RHIC indi- 
cated that the new medium created in ultrarelativistic Au+Au collisions interacted 
even more stron gly than had ronic matter. Moreover, in line with earlier theoretical 
studies in Refs. ^^ '^ ^ the medium showed phenomena of an almost perfect 
liquid of partons l ^ l ^ l, i.e. a strongly coupled system showing high collectivity. The 
most intriguing questions have become: what degrees of freedom are relevant in this 
new sQGP state and how is the transition from the partonic liquid to the gas of 
interacting hadrons realized? 

The c, c quark degrees of freedom are of particular interest in context with the 
phase transition to the sQGP. The heavy flavor sector is important for unraveling 
the actual dynamics from the experimental side since the high masses of the charm 
and especially bottom quarks provide independent (and new) energy scales. The 
hadronic bound states (with a c or c quark) have a much larger mass than the 

that cc pair s can only be formed in the v ery 

25 26 27 28 A s h as hoon nrnnnsprl in Hpf EH 



ordinary hadrons, and it is expected that cc pair s can only be formed in the ver 
early phase of the heavy-ion collision I 25 | 26 | 27 [28| DCCn proposed in Ref. 



the strangeness degrees of freedom might play an important role in distinguishing 
hadronic and partonic dynamics. The kaons (and antikaons), indeed, have been 
proven to provide a suitable probe for the compression phase in heavy-ion reactions 
at SIS energies. One expects that, due to the new scale introduced by the charm 
flavor, Z?-mcsons can be used in a similar way to probe the high density matter 
at FAIR around the threshold for open charm production ^01. Moreover, heavy 
quarkonia (cc, 66) might no longer be bound at high temperature. Due to different 
binding energies, quarkonium states should 'melt' at different temperatur es a nd 
might provide the cleanest "measurement" for the energy densities achieved ^1 

Apart from the total and relative abundances of charmonia and open charm 
mesons also their collective properties are of interest. Here especially the trans- 
verse momentum (or mass) spectra are expected t o Provide v aluable insight into 
the dynamics in cither the very early or late phase | 31 | 32 | 33 |34] ^ e reca n that the 
collective properties of open charm mesons go along with a suppression (relative 
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to scaled pp reactions) of high transverse momen tum px particles as observed at 
RHIC by the STAR and PHENIX Collaborations S3 as we ll as the suppression of 
charmonia (J/v?, VP'). Since the charm quarks are produced early in the reaction, 
their rescattering - reflected in high pt suppression - and collectivity (in the el- 
liptic flow V2(pr)) signal more sensitively the dynamics of the early phase. All the 
arguments above qualify charm-flavored particles as practical and promising probes 
for an exploration of QCD matter. 

The present review is structured as follows: We shall start in Section [2] with a 
description of elementary reactions involving 'charm' and introduce the transport 
models that will be employed for a dynamical description of proton-nucleus and 
nucleus- nucleus collisions, respectively. The latter dynamical models allow for an 
overview of the local energy densities achieved in Au+Au collisions from FAIR to 
RHIC energies (Section[3]). Wc continue with 'cold nuclear matter' effects in proton- 
nucleus reactions in Section |4] and provide information about the 'hadronic envi- 
ronment' of charmonia encountered in relativistic heavy-ion reactions in Section [5l 
The leading concepts/models for charmonium dynamics are presented in Section 6 
as well as their actual implementation in the microscopic transport approach. The 
relative abundances of charmonia and open charm mesons are presented in Sec- 
tion 7 together with a discussion of the anomalous suppression of charmonia ( J/\I/, 
ty') in comparison to the experimental data available. The issue of charmonium 
chemistry is examined in Section [8] in comparison with the statistical hadronization 
model. Sections l9l and [TOl are devoted to the energy loss of charmonia in the dense 
medium, i.e. in particular to the transverse mass spectra and high pt quenching. 
The collective flow properties of open charm and charmonia are studied for the 
elliptic flow vi in Section 1111 We conclude with a summary and discussion of open 
questions in Sections [T2l and [T^l 

2. Transport models and elementary reactions 

The main difficulty in the study of nuclear matter under extreme conditions as 
created in relativistic heavy-ion collisions is that the information about the ini- 
tial sQGP stage of matter can be obtained only indirectly from the measure- 
ment of hadronic observables. The sQGP signal might be strongly distorted by 
the hadronization process and final state interactions of the hadrons. In order to 
reliably subtract the "trivial" hadronic contribution from the sQGP signal, one 
needs a microscopic transport dynamical approach. Also, the possible equilibration 
of quark-gluon matter is on central interest. Such issues of equilibration phenom- 

ena are tradition ally examined within noncquilibrium relativistic transport theory 
Ia6ll7l37l38i:ffll40l 

Ideally, one would use a microscopic transport model containing the proper de- 
grees of freedom - quarks and gluons in the initial phase and hadrons in the final 
phase - and parton dynamics which are consistent with the Lattice QC D equation 
of state. The development of such approaches is in progress | 41 | 42 | 43 | However, a 
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proper understanding of the transport properties of the partonic phase is still lack- 
ing (sec Section [13] for details). One has to keep in mind that most of the results in 
this review have been obtained based on hadron-string or hydrodynamical models 
without including an explicit phase transition to the QGP. As wc will sec, there are a 
number of questions that can be answered before the development of a nonperturba- 
tive parton-hadron transport approach is completed. For example, the interactions 
of J/vfr's with formed hadrons in the late stages of the collision (when the energy 
density falls below a critical value of about 1 GeV/fm 3 corresponding roughly to 
the critical energy density for a parton/hadron phase transition) gives a sizable 
contr ibution to its ano malous suppression at all beam energies as demonstrated in 
Refs. I 44 | 45 | 46 | 33 | 47 | 48 ] Accordingly, this more obvious hadronic contribution has 
to be incorporated when comparing possible models for QGP-induccd charmonium 
suppression to experimental data. 

In any case, the link of the underlying physics to the heavy ion experi- 
ment s shou ld be pr ovided by dynamical trans port m odels, such as UrQMD I 40 | 49 | 
HSD [50136| GiBUU [51152] R q MD [53] q GSM 02I55] or AMPT [42] we recall that 

microscopic (pre-)hadronic transport models describe the formation and distribu- 
tions of many hadronic particles from SIS to SPS energies reasonably well 
Furthermore, the nuclear equation of state has been extracted by comparing tra ns- 
port calculations to hadron flow data up to AGS energies I 57 | 58 | 59 | 60 | 61 | 62 | In 
particular, at SIS energies, microscopic transport models reproduce the data on 
the excitation function of the proton elliptic flow v-i when incorporating a soft, 
momentum-dependent equation of state lfi SI 64 | 65 | 52 [ Tn addition to nucleus-nucleus 
collisions from SIS to SPS energies I 56 | 66 | ; the HSD transport approach is found to 
work reasonably well also at RHIC energies for the 'soft' hadron abundances, so 
that the 'hadronic envi ronment' for open charm mesons and charmonia should be 
sufficiently realistic S3. One also finds generally a good agreement among the pre- 
dictions of different transport models for hadron multiplicities. Modest differences 
between the HSD and UrQMD transport results for pion and kaon multiplicities can 
be attributed to different implementations of string formation and fragmentation, 
which are not sufficiently controlled by experimental data. 

The precision of such models depends crucially on the elementary input, i.e. 
the knowledge of the elementary reaction cross sections. In order to tackle the pro- 
duction and dynamics of charmed hadrons in heavy-ion collisions, one needs data 
from NN and irN reactions at different energies to establish the underlying cross 
sections and distributions. One can apply perturbative QCD (pQCD) to calculate 
the tot al c ross section for cc pair production. The results of next-to-leading order 
pQCD ESI calculations are shown in Fig. Q] in comparison to pp data as well as to 
pA (AA) measurements divided by the number of binary collisions N co u (charm 
pairs, due to their high mass, are expected to be produced only in initial nucleon- 
nucleon collisions and their yield should scale with N co u). The discrepancy of a 
factor of ~ 2 between the PHENIX and STAR cross sections is not yet understood. 
Moreover, the large theoretical uncertainty implies that there is little predictive 
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E769 p+A; PR 433. 127 
NA16 p+A: PR 433, 127 
NA27 p+A: PR 433. 127 
E743 p+A: PR 433. 127 
E653 p+A: PR 433. 127 
HERA-Bp+A: PR 433, 127 
NA50 p+A: PR 433. 127 
NA60 In+ln (preliminary) 



STAR d+Au: PRL 94. 062301 
STAR A+Au (preliminary) 
PHENIX p+p: PRL 97. 252002 
PHENIX Au+Au: PRL 94. 082301 
PHENIX Au+Ali; PRL 88, 192303 




NLOwith CTEQ6M 
NLO Uncertainty Bound 

(from R. Vogt, arXiv:0709.2531) 

i i i i i i i i I i 



10 J 



Vs (GeV) 



Fig. 1. Comparison of total cross section measurements for cc pairs as a function of the invariant 
energy y/s. The STAR and PHENIX results are given as cross section per binary collisions; vertical 
lines reflect the statistical errors, horiz ont al bars indicate the systematic uncertainties (where 
available) . The figure is taken from Ref . ISZJ 



power in the pQCD total charm cross section. Therefore, the use of phenomcno- 
logical parametrizations of the world data for the total charm cross sections from 
elementary reactions is legitimate so far. This approach is followed up in transport 
models. 

In the rest of this Section we present the implementation of charm produc- 
tion from elementary hadron collisions in the Hadron-String-Dynamics (HSD) 
transport modeQ. The total cross sections for the elementary production channels 
including the charmed hadrons D, D, D*, D*, D„ p«, p**, ff* , J "/VP \P(25), Ysv 
from NN and -kN collisions were fitted in Refs. | 36 | 44 | 33 | 47 | 1 8|45| to pyTHIA ESl 
calculations above y/s = 10 GeV and extrapolated to the individual thresholds, 
while the absolute strength of the cross sections was adjusted to the experimental 
data (cf. Ref. ^j). The actual results are displayed in Fig. [5] for p + N and n + N re- 
actions. The total charmonium cross sections (i = \c J/^, from NN collisions 



a The open source code is available from Ref. 
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Fig. 2. The cross section for D + D, J/ty and \E>' meson production in pN (left part) and irN 
reactions (right part). The solid lines show the paramctrizations used in HSD, whereas the symbols 
stand for the experimental data from Refs.l 69 ! 70 | 71|72 | 73|74 | 75 | 76 | 77 | 78 | Note that thc cross 
sections include the decay from \c mesons. The figure is taken from Ref. HH 



as a function of the invariant energy y^s are approximated by the expression 

"!""(>) =fia (1) 

where m, denotes the mass of charmonium i while y/sol = mj+2mjv is the threshold 
in vacuum. The parameters in Eq. (fTJ) have been fixed to describe the J/^> and \P' 
data up to the RHIC energy ^ = 200 GeV (c/. Ref. We use a = 0.16 mb, a 
= 10,j3 = 0.775. 

The parameters fi are fixed as f Xc = 0.636, fj/y = 0.581, = 0.21 in order 
to reproduce the experimental ratio 

B( Xcl -» JMa^ + B {Xc2 -» J/*K c2 = Q 344 ± Q Q31 

measured in pp and ttN reactions I 83 | 84 | as we jj as |- ne av eraged and pA ratio 
(B /J/J (^ r ')cr^/)/(S M/J (J/4 r )<Tj/*) ~ 0.0165 (c/. thc compilation of experimental data 
in Ref. ES) The experimentally measured J/\& cross section includes thc direct J/ty 
component (ctj/*) as well as the decays of higher charmonium states Xc and *S?' , 
i.e. 

"7/1 = <V/« + B (Xc ~> J/*Ke + - W< ■ (2) 

Note, we do not distinguish the Xei(l-P) and x C 2(lP) states. Instead, we use only 
the Xci(l-P) state (which we denote as x c ), however, with an increased branching 
ratio for the decay to J/ty in order to include the contribution of Xc2(l-P), i-e. 
B(Xc -> J/*) = 0.54. Further more, we adopt B{W -» J/*) = 0.557 from Ref.^. 

We recall that (as in Refs. I 33 | 47 | 87 | 88 | 89 |^ charm degrees of freedom in the 
HSD approach are treated perturbatively and that initial hard processes (such as cc 
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pp->J/>P+X, s =200 GeV 




Fig. 3. Cross section for the differential J/* production in rapidity (times the branching ratio to 
di-muons B^^) in pp collisions at y/s = 200 GeV. Th e HS D (input) parametrization (s olid line) 
is compared to the PHENIX data (symbols) from Ref . EH The figure is taken from Ref. EH 

or Drell-Yan production from NN collisions) are 'pre-calculatcd' to achieve a scal- 
ing of the inclusive cross section with the number of projectile and target nucleons 
as Ap x At when integrating over impact parameter. For fixed impact parameter 
b, the cc yield then scales with the number of binary hard collisions N co u ( cf. Fig. 8 

IOOI 

in Ref. To implement this scaling, we separate the production of the hard and 
soft processes: The space-time production vertices of the cc pairs are 'precalculatcd' 
in each transport run by neglecting the soft processes, i.e. the production of light 
quarks and associated mesons, and then reinserted in the dynamical calculation 
at the proper space-time point during the actual calculation that includes all soft 
processes. As shown in Ref. ^Hl this prescription is very well in line with Glauber 
calculations for the production of hard probes at fixed impact parameter, too. We 
mention that this 'precalculation' of cc production might be modified at RHIC en- 
ergies due to changes of the gluon stru cture functions during the heavy-ion reaction 
or related shadowing phenomena ^21. The amount of shadowing at RHIC energies 
will be discussed in more detail in Section |U 

In addition to primary hard NN collisions, the open charm mesons or charmonia 
may also be generated by secondary meson-baryon (mB) reactions. Here we include 
all secondary collisions of mesons wi th b aryons by assuming that the open charm 
cross section (from Section 2 of Ref. only depends on the invariant energy -y/s 
and not on the explicit meson or baryon state. Furthermore, we take into account 
all interactions of 'formed' mesons - after a formation time of Tp w 0.8 fm/c (in 
their rest frame) - with baryons or diquarks. For the total charmonium cross 
sections from meson-baryon (or ttN) reactions we use the parametrization (in line 
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p T 2 [GeV/c] p T 2 [GeV/c] 



Fig. 4. Differential cross section of J/^S production in pp collisions at y/s = 200 GeV at mid- 
rapidity (\y\ < 0.35, l.h.s.) and at forward rapidity (1.2 < \y\ < 2.2, r.h.s.) vs the transverse 
momentum squared p^, as implemented in HSD (solid line) compared to the PHENIX data from 
Ref.El(dots). The figure is taken from Ref.ES 



with Rcf. 

< N (s)=fib (l-^§)' (3) 

with 7 = 7.3 and b = 1.24 mb, which describes the existing experimental data at 
low y/s reasonably well, as seen in Fig. [2] 

Apart from the total cross sections for charmonia we also need the differen- 
tial distributions of the produced mesons in the transverse momentum px and 
the rapidity y (or Feynman xp) from each individual collision. We recall that 
xp = p z /p™ ax ~ 2p z /y/s with p z denoting the longitudinal momentum. For the 
differential distribution in xf from NN and irN collisions we use the ansatz from 
the E672/E706 Collaboration ^3 and for the px distribution: an exponential 

dN 

(1 - \x F \Y exp{-b PT p T ) (4) 



dxpdpj 



with b PT = 2.08 GeV 1 at SPS energies or a power law parametrization from Ref.l^il 
which has been fixed by the STAR data ^] at high energies relevant to RHIC, i.e. 

v (1 _ Mr f 1 + |LV" (5) 



dxFdpr \ b PT / 

with b PT = 3.5 GcV/c and c PT = —8.3. The exponent c is given by c = a/(l + 
b/y/s) and the parameters a, b are chosen as o,nn — 16, bNN = 24 9 G eV for NN 
collisions and a^N = 4.11, b^N = 10.2 GeV for ttN collisions SZE31. Not e that 
the parametrizations of the differential cross sections are taken as in Refs. I 33 | 4T |^ 
apart from a readjustment of the parameter ajvw i n order to reproduce the recently 
measured rap idity distribution of <//<]/ 's in p + p reactions at yfs = 200 GeV by 
PHENIX ^Tin Ref.^Sl. 



August 11, 2008 12:47 WSPC/INSTRUCTION FILE Linnyk 



Charm in heavy ion collisions 11 

The resulting rapidity distribution for J production in pp collisions at ^fs = 
200 GeV is shown in Fig. [3] We also present the pp — > J/^S + X differential cross 
section in p\ at mid-rapidity < 0.35) and at forward rapidity (averaged in the 
interval 1.2 < \y\ < 2.2) in Fig. HI Both kincmatical distributions (in y and pr) are 

Ol I 

in line with the data from Ref. within error bars. 

3. Energy-density evolution in relativistic heavy-ion collisions 

The HSD approach ^ provides the space-time geometry of nucleus- nucleus reac- 
tions and a rather reliable estimate for the local energy densities achieved, since the 
production of secondary particles with ligh t and a single strange quark/ antiquark 
is described well from SIS to RHIC energies | 56 | 66 | ^ gee a j g0 s ec y on 5) j n trans- 
port approach the local energy density is calculated from the energy-momentum 
tensor T Mt/ (x) for all space-time points x in the local rest frame: e(x) = Tqq c (x), 
where Tqq c (x) is calculated from T llv (x) by a Lorentz boost to the local rest frame. 
In order to exclude contributions to from noninteracting nucleons in the initial 
phase all nucleons without prior interactions are discarded in the rapidity intervals 
[ytar - 0.4, y tar + 0.4] and [y pro - 0A,y pro + 0.4] where y tar and y pro denote pro- 
jectile and target rapidity, respectively. Note that the initial rapidity distributions 
of projectile and target nucleons are smeared out due to Fermi motion by about 
Ay w ±0.4. Some comments on the choice of the grid in space-time are in order 
here: In the actual calculation (for Au+Au collisions) the initial grid has a dimen- 
sion of 1 fm x 1 fm x l/7 cm fin, where 7 cm denotes the Lorentz 7-factor in the 
nucleon-nucleon center-of-mass system. After the time of maximum overlap t m of 
the nuclei the grid-size in beam direction Azo = l/7cm [fm] is increased linearly 
in time as Az = Azo + a{t — t m ), where the parameter a is chosen in a way to 
keep the particle number in the local cells of volume AV(t) = AxAyAz(t) roughly 
constant during the longitudinal expansion of the system. In this way local fluctua- 
tions of the energy density e(x) due to fluctuations in the particle number are kept 
low. Furthermore, the time-step is taken as At = 0.2Az(t) and increases in time 
in analogy to Az{t). This choice provides a high resolution in space and time for 
the initial phase and keeps track of the relevant dynamics throughout the entire 
collision history. 

3.1. SPS energies 

As a first example we display in Fig. [5] the energy density s(x,y = 0, z;t) for a 
Pb+Pb collision at 160 A-GeV and impact parameter b = 1 fm in terms of contour 
lines for times of 1, 2, 3 and 5 fin/c (from contact). It is clearly seen that energy 
densities above 4 GeV/fm 3 are reached in the early overlap phase of the reaction 
and that e(x) drops within a few fm/c below 1 GeV/fm 3 in the center of the grid. On 
the other hand the energy density in the region of the leading particles - moving 
almost with the velocity of light - stays above 1 GeV/fm 3 due to Lorentz time 
dilatation since the time t here is measured in the nucleon-nucleon ccntcr-of-mass 
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Fig. 5. The energy density e(x,y = 0, z;t) from HSD for a Pb+Pb collision at 160 A-GeV and 
impact parameter b = 1 fm in terms of contour lines (0.01, 1, 2, 3, 4 GeV/fm 3 ) for times of 1, 2, 3 
and 5 fm/c (from contact). Note that noninteracting nuclcons have been discarded in the actual 
calculation of the energy-momentum tensor. The figure is taken from Ref. H2. 

system. Note that in the local rest frame of the leading particles the eigentime r is 
roughly given by r m t/-f cm with j cm w 9.3 (at 160 A-GeV). 

Another view of the space time evolution of the energy density is given in Fig. [6] 
where we display e(x = 0, y = 0, z; t) for the same system as in Fig. [5] on a linear 
scale. The contact time of the two Pb nuclei here is 2 fm/c and the overlap phase 
of the Lorentz contracted nuclei is identified by a sharp peak in space-time which 
is essentially given by the diameter of the nuclei divided by 7 cm . As noted before, 
the energy density in the center of the reaction volume (z « 0) drops fast below 1 
GeV/fm 3 whereas the ridges close to the light-cone basically stem from the leading 
ends of the strings formed in the early nucleon-nucleon collisions. In these space- 
time regions all reaction rates are reduced by the factor ~ l/7 cm such that the 
transport calculations have to be carried out to large times of several hundred fm/c 
in order to catch the dynamics and decays in these regions. In the central regime, 
however, all interaction rates vanish after about 15 fm/c. Since the c, c pairs are 
produced dominantly at midrapidity with a small spread in rapidity (<r y m 0.8 at 
160 A- GeV) it is the central region that is of primary interest for this study. 
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Fig. 6. The energy density e(x = 0, y = 0, z;t) from HSD for a Pb+Pb collision at 160 A-GeV 
and impact parameter 6=1 fm on a linear scale. Note that noninteracting nucleons have been 
discarded in the actual calculation of the energy-momentum tensor such that e(x ) ^ only after 
contact of the two Pb nuclei which is ~ 2 fm/c. The figure is taken from Ref. 



3.2. RHIC energies 

The energy density e(r;£) becomes very high in a central Au+Au collision at 
y/s = 200 GcV as shown in Fig. UJ (in analogy to Fig. [6] for the top SPS energies). 
Fig. [7J shows the space-time evolution of the energy density e(x = Q,y — 0, z;t) 
for a Au+Au collision at 21300 AGcV or y/s — 200 GeV. It is clearly seen that 
energy densities above 16 GeV/fm 3 are reached in the early overlap phase of the 
reaction and that e(x) drops after about 6 fm/c (starting from contact) below 1 
GeV/fm 3 in the center of the grid. On the other hand the energy density in the 
region of the leading particles - moving almost with the velocity of light - stays 
above 1 GeV/fm 3 due to Lorentz time dilatation since the time t in the transport 
calculation is measured in the nucleon-nucleon center-of-mass system. As seen from 
Fig. [Jj the energy density in the local rest frame is a rapidly changing function 
of time in nucleus- nucleus collisions. For orientation let us quote the relevant time 
scales (in the cms reference frame): 

- The cc formation time r c w 1/M± is about 0.05 fm/c for a transverse mass 
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Fig. 7. The energy density e(x = 0, y = 0, z; from HSD for a central Au+Au collision at y/s = 
200 GeV. The time t is given in the nucleon-nucleon center-of-mass system. The figure is taken 
from Ref . 



of 4 GeV; the transient time for a central Au+Au collision at y/s = 200 GeV is 
t r w 2i? J 4/7 cm « 0.13 fm/c. According to standard assumptions, the cc pairs are 
produced in the initial hard NN collisions dominantly by gluon fusion in the time 
period t r . In fact, the formation time r c is significantly smaller than t r , which 
implies that the c or c quarks may interact with the impinging nucleons of the 
projectile or target for times t < t r . 

- Using the Bjorken estimate for the energy density and employing the time- 
scale t r = 0.13 fm/c, the energy density - after the nuclei have punched through 
each other - amount s to ab out 5/0.13 > 30 GeV/fm 3 (as quoted also in the HSD 
calculations in Refs. | 45 | 82 fy g; ven w hen adding the cc formation time, this gives 
an energy density ~ 5/0.18 « 28 GeV/fm 3 . So the numbers in Fig. [7] agree with 
transparent and simple estimates and illustrate the high initial densities after cc 
production from primary interactions. 

The energy densities quoted above are considerably different from the estimate 

<E T > f- 

T ■ ^ = (6) 

where < Et > is the average transverse energy per particle, dN / drj the number of 
particles per unit of pscudorapidity, and r a formation time parameter often used 
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Fig. 8. Left part: The transverse energy Et per pseudorapidity interval dn divided by the nu mber 
of participant pairs (0.5N par t) from HSD (solid line) in comparison to the PHENIX data (dots)EEl 
Right part: The Bjorken energy density esj • t from HSD (solid line) for Au+Au collisions at y/s 
= 200 GeV in comparison to the PHENIX data (dots) ED The figure is taken from Ref . IS21. 



as r = 1 fm/c. Furthermore, irR^ denotes the overlap area for the corresponding 
centrality. Is is important to point out that the estimate ^ is only well defined 
for the product t • esj! The question naturally arises, if the transport calculations 
follow the corresponding experimental constraints. 



To this aim we show dEj< /drj (divided by half the number of participants N p 
from HSD in Fig. [8] (l.h.s.) in comparison to the measurements by PHENIX 




Accordingly, the Bjorken energy density tBj - multiplied by the time-scale t 
from HSD is shown additionally in the r.h.s. of Fig.[5]in comparison to the PHENIX 
data as a function of N part . The similarity between the calculated quantities and 
the experimental data demonstrates that the space-time evolution of the energy- 
momentum tensor T^ v in HSD is sufficiently well under control also at RHIC ener- 
gies. 

3.3. FAIR energies 

The question emerges if central collisions of e.g. Au+Au at the future FAIR facil- 
ity might be also suited to explore a possible phase transition to the sQGP. For 
a quantitative orientation we display in Fig. [5] the volume (in the nucleus-nucleus 
center-of-mass) with energy densities above 1 GeV/fm 3 and 2 GeV/fm 3 as a func- 
tion of time for a central Au+Au collision at 25 A-GcV, where only interacting and 
produced hadrons have been counted as explained in the beginning of this Section. 
It is important to note that in HSD the high energy density is essentially build up 
from 'strings', i.e. 'unformed' hadrons. The absolute numbers in Fig. [9] have to be 
compared to the volume of a Au-nuclcus in the moving frame which - for a Lorentz 
7-factor of 3.78 - gives f» 330 fm 3 . In this case, the overlap phase of the nuclei 
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Fig. 9. Time evolution for the volume with energy density e > 1 GeV/fm 3 (dashed line) and > 
2 GeV/fm 3 (solid line) i n th e HSD approach for a central Au + Au reaction at 25 A-GeV. The 
figure is taken from Ref. 1221. 

lasts for about 3.7 fm/c during which energy densities above 2 GeV/fm 3 are seen 
in Fig. [9] in a sizeable volume. Thus also at 25 A-GeV the phase boundary to a 
QGP might be probed in a sizeable volume for time scales of a few fm/c. Contrary 
to central collisions at the SPS and RHIC, these volumes are characterized by a 
high net quark density; for such configurations we presently have only 'hints' from 
lattice QCD calculations rather than solid information. 

4. Proton-nucleus reactions: Cold nuclear matter effects 

Before coming to charm production and propagation in heavy-ion reactions it is 
mandatory to explore the charm dynamics in proton-nucleus reactions. Such reac- 
tions are a first step beyond the elementary pp reactions and provide an additional 
reference with respect to the heavy-ion case. Since inp+A reactions only subnormal 
nuclear densities are achieved and the target nucleus remains approximately in its 
geometrical shape for the first few fm/c, these reactions are much easier to treat and 
allow to separate 'cold nuclear matter' (CNM) effects from those induced by the 
new partonic medium encountered in relativistic nucleus-nucleus collisions. As we 
will show in the following, one can observe e.g. gluon shadowing at RHIC energies 
by it's influence on charmonium production in d + Au collisions. Additionally, the 
amount of 'normal' nuclear charmonium dissociation by the target nucleons can be 
probed. Thus p + A reactions provide a necessary base-line for the heavy-ion case 
and can independently be controlled by experimental data. 

It is found experimentally that the yield of J/^ in p + A and A + A reactions 
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Fig. 10. The ratio <r j /^/cdy (times the branching ratio -B^) measured in p-A and nucleus- nucleus 
collisions at the SPS, compiled and re-scaled, when necessary, to 158 GeV incident energy by the 
NA60 Collaboration. The quantity L denotes the effective path length of the J/\E> as extracted 
within the Glauber model. The lines indicate the results of a Glauber fit to the p-A data and the 
size of the error as calculated by the NA60 Collaboration. The figure is taken from Ref. HI. 

is modified compared t o tha t in p + p collisions - scaled with the number of initial 
binary scatterings iV co ;jlSEZI Indeed, the produced cc can be dissociated or absorbed 
on either the residual nucleons of the projectile or target or on light co-moving 
particles (usually on mesons or, at high energy, on partons) produced in the very 
early phase. Alternatively, the initial production of cc pairs by gluon fusion might 
be suppressed due to shadowing (at RHIC energies). In particular, charmonium 
absorption on baryons is the leading suppression mechanism in d + A (p + A) 
scattering at SPS energies and is an important base-line for the investigation of 
charmonium absorption. 

4.1. SPS energies 

The NA50 and NA60 Collaborations present their results on J/ty suppression as 
the ratio of the dimuon decay of the J/^S relative to the Drcll-Yan background 
from 2.9 - 4.5 GeV invariant mass as a function of the transverse energy Ex, or 
alternative, as a function of the number of participants N pax t, i.e. 

B MM a(J/tf)/a(DY)|a. 8 -4.6, (7) 
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where is the branching ratio for J/^ — > /i + u~ . In order to compare our 

calculated results to experimental data (see below) we need an extra input, 
i.e. the normalization factor B^aNNiJ/ty/PNNiDY), which defines the J/^S 
over Drell-Yan ratio for elementary nuclcon-nucleon collisions. We will adopt 
B^<T NN {J/^)/a NN {DY) = 36 in line with the NA60 compilation El (at the SPS 
energy of 158 GeV). 

In Fig. [10] we show the combined p + A, S + U, Pb + Pb and In + In data (from 
NA50 and NA60) for the ratio B^ L a(J '/*) / 'a(DY) at 158 A- GeV as a function of 
centrality - reflected in the effective path length L - together with the Glauber-model 
fit by the NA50/60 collaboration. The 'default' interpretation of the experimental 
results in Fig. \W\ is that for a J/vP path length L below about 7 fm dominantly 
'normal' J/^S dissociation with target nuclcons is seen while for L > 7 fm an 
'anomalous' suppression sets in. Since this 'anomalous' suppression only is observed 
in central In + In and Pb + Pb collisions it is attributed to a 'hot matter' effect 
in contrast to the 'normal' absorption ('cold nuclear matter' effect). However, in 
order to distinguish more clearly such 'hot' and 'cold nuclear' matter effects it is 
mandatory to employ non-equilibrium transport. 

In order to study the effect of charmonium rescattering on projectile/target 
nuclcons, we adopt in HSD the following dissociation cross sections of charmonia 
with baryons indep ende nt of the energy (in line with the most recent NA50 and 
NA60 compilations 

OccB = 4.18 mb; (8) 

a J/^B = 4.18 mb; o Xc b = 4.18 mb; ct*'_b = 7.6 mb. 

The applicability of the Glauber picture to the baryon-induced suppression (at 
L < 7 fm) as illustrated in Fig. [10] suggests that the produced cc pair can be 
absorbed on baryons already in its pre-rcsonant state. In |8j the cross section a cS B 
stands for a (color dipole) pre-resonance (cc) - baryon cross section, since the cc 
pair produced initially cannot be identified with a particular charmonium due to 
the uncertainty relation in energy and time. For the life-time of the pre-reson ance 
cc pair (in it's rest frame) a value of t cE — 0.3 fm/c is assumed following Ref. 1^21. 
This time scale corresponds to the mass difference of the and J/ty according to 
the uncertainty relation. Note that - in contrast to the absorption on primordial 
baryons B (nucleons of the incoming nuclei) - interactions with secondary particles 
created in the nucleus-nucleus collision (mesons or secondary baryons) are only 
allowed after the local energy-density has dropped below 1 GeV/fm 3 in order to 
assure that the interaction is hadronic. 

In Fig. [TTJ we present the results of the HSD calculations for the observable 
([7]) for Pb+Pb and In+In collisions in the nuclear suppression scenario, i.e. with 
only baryonic absorption and no additional, meson- or parton-induced suppression. 
Instead of the (model-dependent) path length L we display this ratio as a function 
of the number of participants N part , which can directly be taken from the transport 
calculations. The dashed (blue) lines stand for the HSD results while the (green- 
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Fig. 11. The ratio B ll/J ,o-(J/fy)/o-(DY) as a function of the number of participants in In+In (l.h.s.) 
and Pb+Pb reactions (r.h.s.) at 158 A ■GeV. The full symbols denote the data from the NA50 
and NA60 Collaborations (from Rcfs. , while the dashed (blue) lines represent the HSD 

calculations including only dissociation channels with nucleons. The (dashed blue-green) bands in 
the upper parts of the figure give the estimate for the normal nuclear J/ absorption as calculated 
by the NA60 Collaboration. The vertical lines on the graphs reflect the the oretical uncertainty 
due to limited statistics of the calculations. The figure is taken from Ref. E2. 



blue) bands give the estimate for the normal nuclear J/\& absorption as calculated 
by the NA60 Collaboration in the Glauber model^. The normal nuclear suppression 
from HSD is seen to be slightly lower than the (model dependent) estimate from 
NA60, however, agrees quite well with their calculations for more central r eactions. 
The various experimental data points have been taken from Refs. | 4 | 5 | 100 | 

Apart from the statistical uncertainties in the calculations - reflected by the 
vertical lines on the theoretical graphs in Fig. [IT] - some dependence on the model 
parameters enters the actual numbers in Fig. 1111 The charmonium nuclear absorp- 
tion cross section is considered to be 'fixed' by the NA50/NA60 compilations and 
we have taken the same cross section for the 'pre-resonance' cross section for the 
J/^! and Xc- Accordingly, the life-time of the pre-resonance state (r ce = 0.3 fm/c) 
has no impact on the absorption with baryons as far as the J/vI/ and Xc mesons 
are concerned. Only for vp' collisions with baryons this plays a role, since the \&' 
+ baryon cross section is larger (7.6 mb). Consequently, the J/^f suppression (in- 
cluding the feed down from Xc) does not depend on r cg . Within these systematic 
uncertainties we will now be able to separate 'cold nuclear' and 'hot' matter effects 
in relativistic nucleus-nucleus collisions at SPS energies. 
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Fig. 12. Ratio of gluon densities in lead nucleus and in proton as a function of Bjorkcn x at a fixed 
hard scale Q 2 = 5 GeV 2 from various parton distribution fits. The figure is taken from Ref. EH] 



4.2. RHIC energies 

The cross sections <Jj/^Ni&ccN at RHIC energies are currently debated in the lit- 
erature. On one side, all the data on J/ty production in p + A reactions at energies 
y/s < 40 GeV are fo und to be consiste nt with an energy-independent cross section of 
the order of 4 — 7 mbl 5 | 98 | 102 | 103 | 1 04|. Qn other hand, at the much higher energy 
of \fs = 200 GeV some part of the suppression might be attri buted to oth er (initial- 
state) 'cold nuclear matter' effects such as gluon shadowing 1 1 05 | 1 06 | 1 07 | radiative 
gluon energy loss in the initial state or multiple gluon rescattering | 108 | 109 |HQ| 
We recall that 'shadowing' is a depletion of the low-momentum parton distribu- 
tion in a nuclcon embedded in a nucleus compared to the population in a free 
nuclcon; this also leads to a lowering in the (scaled) charmonium production cross 
section in p + A relative to pp reactions. The reasons for depletion, though, are 
numerous, and models of shadowing vary accordingly. There is, therefore, a con- 
siderable (abou^aJ^cjOTof3) uncertainty in the amount of shadowing predicted 
at RHIC I 105 | 106 | 107 | lll | 112 | i nc jeed, there is currently a lack of precise differen- 
tial data, which would allow to constrain the gluon distribution in nuclei at low 
Bjorken x; this region is probed by charmonium production by gluon fusion pro- 
cesses. Fig. [12] shows an overview of the various model concepts which demonstrates 
that at low Bjorkcn x (of the order of 10 _2 -10 -4 ) the nuclear gluon distribution 
function has large uncertainties. 

In the analysis of the d + Au data at y/s = 200 Ge V, in which the maximum 
estimate for the effect of shadowing was made additional absorption 
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Fig. 13. The J/ip nuclear modification factor in d + Au collisions at y/s = 200 GeV as a function of 
rapidity LLUil compared with calculations including shadowing on top of the Glauber model nuclear 
absorption with an adjusted J/tpN breakup cross se ction. Le ft panel: the data are compared with 
calculations HQZl using the EKS98 param eterization | H4 | 11 5] of shadowing. Middle pa nel: sa me as 
the left panel, but the calcu lations Hfi^l employ the nDSg shado wing parameteriza tion ESI. Right 
panel: the calcul ations ll lZl employ the CF picture of shadowing | 107 | 109 | 1 10 |108| r £\ le figures are 
taken from Refs. I 67 l 117 [ 



on baryons allowed by the data was found to lead to uj/^m = 1 — 3 mb or higher, if 
some contribution of anti-shadowing is present. The authors of Ref. advocate 
at least <Jj/^n = 3 mb in order to preserve the agreement with the data of the 
Fermilab experiment E866. The PHENIX Collaboration HEI finds a breakup cross 
section of 2.8^4 mb (using EKS shadowing) which still overlaps with the CERN- 
SPS value of 4.18 mb (though with large error bars). 

Fig. 1131 demonstrates the variation in the rapidity dependence and the amount 
of initial state interaction between different implementations of shadowing for the 
nuclear modification factor defined as 

_ dNfft/dy 
UdA ~ (N coll )-dN^/dy- W 

In Eq. © dNjfig/dy is the J/ty invariant yield in d + A collisions, dN^^/dy is 
the J invariant yield in p + p collisions; (N co u) is the average number of binary 
collisions for the same rapidity bin ((N co u) — 7.6 ± 0.3 according to the PHENIX 
estimate ^-^ j). There is an additional large theoretical uncertainty in the results 
shown in Fig. [13] since in the works above only an approximate model for baryonic 
absorption was applied and not a microscopic transport approach that e.g. also 
includes secondary production channels of charm pairs as described in Section [2] 

We continue our investigation of 'cold nuclear matter' effects at RHIC energies 
for d + Au reactions employing the same cross sections for baryonic absorption ([8} 
as at SPS energies (cf. Fig. [TT]) . In Fig. [14] we compare the HSD result (neglecting 
shadowing) fo r the J/\V production in d + Au collisions at y/s = 200 GeV to the 
PHENIX dataEHI. 

It is seen from Fig.[T5]that the calculations follow approximately 
the decrease in RdA with rapidity, however, with a tendency to overshoot at forward 
rapidity. Within error bars wc find the values of a C cB from ijHJ) to be compatible with 
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Fig. 14. J/ty production cross section in d + Au collisions relative to that in p + p collisions (see 
text for the definit ion o f RdA) tn HSD (red stars) — neglecting gluon shadowing - as compared to 
the PHENIX dataHEH (full dots). The figure is taken from Ref. E21 



the inclusive RHIC measurement as well as with the lower energy data Gil. This 
finding is also in line with the analysis of the PHENIX Collaboration in Ref. HEH. 

In order to shed some further light on the role of shadowing, we compare our 
calculations for RdA in different rapidity bins as a function of the centrality of the 
d + Au collision, which in Fig. [15] is represented by the number of binary collisions 
N co u. The latter number is directly taken from the number of binary hard NN 
collisions in the transport calculation while the comparison with experime nt is 
based on a Glauber model analysis of the data similar to that performed in Ref.^Sl. 
The actual results displa yed in Fig. [15] (stars connected by dashed lines) and the 
PHENIX data from Ref. ^321 are roughly compatible for the rapidity intervals - 
2.2 < y < -1.2 and \y\ < 0.35, but demonstrate that the suppression at forward 
rapidity (1.2 < y <2.2) is underestimated in the color-dipolc dissociation model 
with a constant cross section of 4.18 mb. This clearly points to the presence of 
shadowing effects at least at forward rapidities which is not so pronounced in the 
inclusive data set in Fig. [M] A more serious question is a quantification of the 
shadowing, which is extremely challenging because of the limited statistics of both 
the experimental data and the calculations. Here we do not attempt to attribute 
a fixed number for the shadowing effect but merely point out that independent 
high statistics data will be necessary to fix this unsatisfactory situation from the 
experimental side. 

Nevertheless, some note of caution is appropriate for the further analysis of char- 
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Fig. 15. The ratio R^a ® f° r backward, central and forward rapidity bins as a function of the 
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without including low-x gluon shadowing and slightly overestimate RdA m the forward interval 
1.2 < y < 2.2. The theoretical error bars are due to the finite statistics of the calculation. The 
figure is taken from Ref. 1221 



monium suppression in Au + Au collisions: There are 'cold nuclear matter effects' 
such as 'gluon shadowing' beyond those incorporated in the HSD transport cal- 
culations, and especially quantitative statements about any 'agreement with data' 
might have to be reconsidered. In case of Au + Au reactions the shadowing from 
projectile/target will show up symmetrically around y = and in part contribute 
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to the stronger J/\I/ suppress ion a t forward/backward rapidities. Nevertheless, fol- 
lowing Granier de Cassagnac an anomalous suppression of J/^ beyond 'cold 
nuclear matter' effects is clearly present in the Au + Au data (to be investigated 
below) . 

4.3. FAIR energies 

At FAIR energies shadowing is not expected to show any sizeable effect (as at SPS 
energies). Accordingly, the charmonium suppression will be driven by dissociation 
reactions with baryons, mesons etc. as will be discussed in Section 6. As we will see 
in Subsection 17.51 HSD predicts that the dissociation on nucleons will actually be 
the leading mechanism for J/^i> suppression at FAIR energies. 

5. Hadron abundances from heavy-ion collisions 

We here recall the information gained from previous experimental (and theoretical) 
studies on hadron production in heavy-ion reactions. A general overview on the 
experimental meson and strange baryon abundances from central nucleus-nucleus 
collisions (Au+Au or Pb+Pb) is given in Fig. Ql)l which shows the meson abun- 
dances from central Au+Au reactions as predicted by HSD transport calculations 
in 2000 1331 from SIS to RHIC energies. The experimental data at AGS (dots), SPS 
(squares) and RHIC energies (triangles) have been added recently. All meson mul- 
tiplicities show a monotonic increase with bombarding energy which is only very 
steep at 'subthreshold' energies, i.e. at bombarding energies per nucleon below the 
threshold in free space for NN collisions. Note that the HSD transport results are 
also in a fair agreement with results from the UrQMD transport appr oac h and 
data (from AGS to RHIC energies) for the hadron rapidity distributions ^H. Only 
in case of transverse mass spectra both transport ap pro aches underestimate the 
experimental slopes from lower SPS to RHIC energies^! 

Some comments with respect to the charm production in heavy-ion reactions 
are in order: The mass of the J - for the states of interest here - gives the lowest 
scale of 3.097 GeV. Accordingly, the formation of a J from an initial cc pair 
is the only allowed process (in vacuum) close to the charm threshold, because the 
D + D channel, i.e. N+N^D + D+N + N, has an effective invariant mass of 
3.739 GeV. The associated production of a D(c) meson with a charmed hyperon 
A C (E C ), i.e. N + N — > A c + D + N , is more favorable due to effective invariant mass 
of 3.216 GeV (3.386 GeV), which is lower than for the production of a DD pair. 
This explains, why in Fig. [16] the -D(c) cross section is larger than the D(c) cross 
section close to threshold energies, and the J formation dominates in the far 
subthreshold domain. At higher bombarding energies the meson abundances group 
according to their quark content, i.e. the multiplicities are reduced (relative to ir + ) 
by about a factor of 5 for a strange quark, a factor of r* 2 ■ 5 2 = 50 for ss = <fi, a 
factor of w 100 for D (D), and w 2 • 10 4 for cc = J/*. 
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Energy [A GeV] 

Fig. 16. Overview on the experimental meson abundances and HSD predictions for the multiplic- 
ities of 7r + , n, K + , K~ , <j>, D,D and J/^-mesons for central collisions of Au+Au as a function 
of bombarding energy from SIS to RHIC energies. The figure is taken from Ref. while the 
experimental data have been added recently. 

According to the arguments given above the Z?-mcsons with a c are produced 
more frequently at low energies (due to an associated baryon A C ,E C ). At roughly 
15 A-GeV the cross sections for open charm and charmonia are expected to be of 
similar magnitude, while at higher energies the ratio of open charm to charmonium 
bound states increases rapidly with energy. Since the excitation function for open 
charm drops very fast with decreasing bombarding energy, experiments around 25 
A-GcV, e.g. at the future FAIR facility^, will be a challenging task, because the 
multiplicity of the other mesons is higher by ord ers of magnitude. 

In the statistical hadronization model of Ref. ^2Qi production of cc pairs is 
assumed to proceed by hard initial nucleon-nucleon scattering - as in HSD - but 
the redistribution of charm quarks and antiquarks in the hadronization process is 
assumed to follow statistical laws determined by a chemical freeze-out tempera- 
ture T c f r and the baryon chemical potential [is- The latter parameters are taken 
from th e exp erimental systematics of the chemical freeze-out line in the (T,/ib) 
diagram EDI. The resulting energy dependence of charmed hadron production in 
heavy-ion collisions (at midrapidity) is displayed in Fig. I17| where the absolute 
yields arc shown on the l.h.s. whereas the yields relative to the number of cc pairs 
produced are given in the right panel. Note (in both panels) the scale factors of 
10 and 100 for J/ty and 'J' mesons, respectively. The relative ratios of charmed 
baryons to charmed mesons differ significantly from those in the HSD approach 
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Fig. 17. Energy dependence of charmed hadr qn pr oduction in heavy-ion collisions at midrapidity 



in the statistical hadronization model of Ref. ' u l. Left panel: absolute yields, right panel: yields 
relative to the number of cc pairs. Note (in both panels) t he sc ale factors of 10 and 100 for J/ty 
and mesons, respectively. The figure is taken from Ref. 1121]. 



where full chemical equilibrium is not achieved in the charm sector (cf. Section 8). 
On the other hand it is unknown whether the statistical model is applicable to 
the charm sector in the FAIR energy range (particularly for non-central collisions). 
Future experiments at FA IR a re expected to clarify this issue. 

As pointed out in Ref. E21 dropping D, D masses with baryon density (and/or 
temperature) might lead to an increase of J absorption and to a net lowering of 
the ty'/J/ty ratio for central collisions. Thus the , ^ ^/ /J/\I , ratio could also qualify as 
a probe of _D-mcson in-medium effects. On the other hand, medium modification 
of D-mesons is more complex than a simple picture of dropping masses suggests. 
Microscopic G-matrix studies indicate that the drop of the D-meson masses with 
baryon density is o nly very moderate and as a leading effect one should expect a 
spectral broadening I 123 | 1 24] 

Note, however, that the elementary cross sections for open charm and charmonia 
in pp and nN reactions have to be measured in the relevant kinematical regimes 
before reliable conclusions can be drawn about charm dynamics in the nucleus- 
nucleus case. Experimental data in the 20 - 30 A- GeV with light and heavy systems 
will have to clarify, furthermore, if the quasi-particlc picture of open charm mesons 
at high baryon density is applicable at all or if the dynamics is already governed 
by partonic degrees of freedom rather than hadronic ones (see below). 

The production of 'ordinary' hadrons (with 'u,d,s' flavor) as well as open 
charm and char moniu m at SPS and RHIC ener gies has been calculated so far 
within the AMPT^l, HSD ^^^^ and UrQMD^ll 

transport approaches using 
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paramctrizations for the elementary production channels as described in Section [2j 
Backward channels 'charm + anticharm meson — > charmonia + meson' are treated 
in HSD via detailed balance in a schematic interaction model with a single matrix 
element | M| 2 that is fixed by the J/_$f suppression data from the NA50 collabora- 
tion G^S a t SPS energies (cf. Ref. 122 and Section 6.1). The independent transport 
approaches provide in general very similar results for energy densities, baryon densi- 
ties, meson densities etc. such that the bulk dynamics of relativistic nucleus-nucleus 
collisions is known to a sufficient extent (cf. Ref. for a more detailed comparison). 



6. Anomalous suppression of J/*Sf: The basic models 

In the past, the charmonia J/^ff, Xc, W have been discussed in context of the phase 
transition to the QGP, since cc states might no longer be formed due to color 
screening 1^11. However, more recent calculations within lattice QCD (1QCD) have 
shown that at least the J/"ff survives up to ~ 1.5 T c (T c f=a 0.17 - 0.19 GeV) such 
that the l owest cc states remain bound up to energy densities of about 5 GeV/fm 3 
(see Refs. H2. 9 | 130 fy j s presently not clear, if also the D or D* mesons survive at 
temperatures above T c , but strong correlations between a light quark (antiquark) 
and a charm antiquark (quark) are likely to persist. One may speculate that similar 
correlations survive also in the light quark sector above T c , such that 'pre-hadronic 
comovers' - most likely with different spectral functions - might show up also at 
energy densities above 1 GeV/fm 3 , which is taken as a characteristic scale for the 
critical energy density. 

On the other hand, it is well known that the baryonic (normal) absorption alone 
cannot ex plain the suppression of charmonia in heavy-ion collisions with increasing 
centrality^^ (cf. Fig. fTTj) . Different mechanism for the additional (anomalous) sup- 
pression or formation of charmonia have been suggested in the pa st, i .e. charmonia 
might be 'melting' according to the scenari o ad vocated in Ref. I^SI (their forma- 
tion be suppressed due to plasm a screening ^ absorbed on co- moving mesons 
in hot hadronic matter I H2 | 133 | Qr they could be absorbed early by neighboring 
strings Moreover, charmonia might also be generated in a statistical fashion 
at the phase boundary between the QGP and an interacting hadron gas such that 
their abundance could be in statistica l (chemic al) equilibrium with the light and 
strange hadrons as suggested in Refs. | 134 | 135 | j a tter picture is expected to 
lead not to a suppression but to an enhancement of J /iff mesons at the f ull RHIC 
energy if compared to the scaled J /iff multiplicity from pp collisions 



6.1. 'Comover' suppression (and recombination) 

First of all let us stress that the interactions with 'comoving' mesons lead not only 
to the dissociation of charmonia, but also to their recreation via the inverse recombi- 
nation process D + D — > cc + m, where m — {ir, p, u>, K, ...}. As already pointed out 
before, the J/^f, Xc, formation cross sections by open charm mesons or the inverse 
'comover' dissociation cross sections are not well known and the significance of these 
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channels is discussed controversially in the lite rature | 132 | 137 | 136 | 138 | 139 | 140 | 141 | 
In HSD - following the concept of Refs. | 47 | 48 | _ a s i m pl e 2-body transition model 
is employed with a single parameter \Mq\ 2 that allows to implement the backward 
reactions uniquely by exploiting detailed balance for each individual channel. We 
briefly review this concept in the following. 

Since the charmonium-meson dissociation and backward reactions typically oc- 
cur with low relative momenta ('comovers'), it is legitimate to write the cross section 
for the process 1 + 2^3 + 4 as 

cr 1+W4 ( S ) = 2^^^±\M^ ( T^iVEL, (10) 

where Ek denotes the energy of hadron k (k = 1,2,3,4), respectively. The initial 
and final momenta for fixed invariant energy y/s are given by 



2 (s - (mi + m 2 Y)(s - (mi - m 2 ) ) 
Pt = 



4s 

2 



_ (s - (m 3 + m 4 ) 2 )(s - (m 3 - to 4 ) 2 ) 



where m,k denotes the mass of hadron fc. In (|10[) \Mi\ 2 (i = Xc, J/ ^ stands 
for the effective matrix element squared, which for the different 2-body channels is 
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Fig. 19. The distribution in the invariant collision energy y/s for J/^S absorption on 7T-, 1J-, p- 
and u-mesons in central (b = 2 fm) Pb + Pb collisions at 160 A- GeV within the HSD transport 
approach. The figure is taken from Ref. I2£J 



taken of the form 

\M l \ 2 = \M t \ 2 for {■K,p) + {cc) i ^D + D (12) 
\Mi\ 2 = 3\Mi\ 2 for (ir,p) + {cc) l -► D* +D, D + D* , D* + D* 

\M t \ 2 = i|A/ ? | 2 for (A', AT*) + (cc)< As +5, 5 S + D 

\M t \ 2 = |AA| 2 for (A, A*) + (cc)< -» D s + 5*, 5, + D*, D* s + D, 

D* s +D, D* + D* 

The relative factors of 3 in (fT2")) arc guided by the sum rule studies in Ref. 11421 
which suggest that the cross section is increased whenever a vector meson D* or 
D* appears in the final channel while another factor of 1/3 is introduced for each s or 
s quark involved. The factor ((7713 + TO4)/ v / s) 6 in (fT0|) accounts for the suppression 
of binary channels with increasing yfs and has been fitted to the expe rimental data 
for the reactions 7r + N -> p + N, u + N, cj) + N, K + + A in Ref. 1421. 

We use the same matrix elements for the dissociation of all charmonium states 
i (i = Xc, J with mesons: 

\Mj/v\ 2 = \AI x f = \M*,\ 2 = \M \ 2 . (13) 

The best fit for |M | 2 (in comparison to the latest NA50 and NA60 analysis EG^) 
is obtained for |A/ | 2 = 0.18 fm 2 /GcV 2 SI this value will be employed in the HSD 
calculations for all bombarding energies and systems. 

The resulting J /iff dissociation cross sections with ir, p, K and A* mesons 
are displayed in Fig. [18] . The final state includes all binary channels compatible 
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Fig. 20. The cross sections for the channels D+D, D+D* , D*+D, D*+D* — > J/* + meson (upper 
part) and the channels involving s or s quarks D s +D, D s +D*, D*+D, D*+D+ — > J/ty + (K, K*) 
(lower part) as a function of the invariant energy yfs according to the model described in the text 
(Eq. Gil). 



with charm quark and charge conservation. Note that for the comover absorption 
scenario essentially the regime 3.8 GcV < ^fs < 4.8 GcV is of relevance (cf. Fig.ll9|) 
where the dissociation cross sections are on the level of a few mb. We note that 
the expli cit channel J/ 1 ^ + it — > D + D, which has often been calculated in the 
literature ^^^^^^^^^^^ , is below 0.7 mb in our model. A somewhat more essential 
result is that the J/ty dissociation cros s section with p-mesons is in the order of 5-7 
mb as in the calculations of Hag linESl used before in Ref. since this channel was 
found to dominate the J/fy dissociation at SPS energies^^Sl. Indeed, the calculations 
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Fig. 21. The J/ty dissociation cross section specified in 11101 -1 1121 1. uppe r blue line, compared to 
the exclusive cross sections with DD*, D* D* final states from Ref. | 132 | ( rec [ anc j g reen lines). 



show that in the realistic fireball evolution the J/'t's are dominantly absorbed in 
collisions with p-mcsons as seen in Fig. 1191 where the distribution in the invariant 
collision energy y/s is plotted for J/^S absorption on 7T, 77, p and oj mesons in central 
(6 = 2 fm) Pb+Pb collisions at 160 A-GeV. It should be pointed out that the 
'comovcr' dissociation channels for charmonia are described in HSD with the proper 
individual thresholds for e ach chan nel in contrast to the more schematic 'comover' 
absorption model in Refs. | H2 | 133 |_ 

The advantage of the model (jlOjl is that detailed balance for the binary reactions 
can be employed strictly for each individual channel, i.e. 



| (2S 1 + 1)(2S 2 + 1) pf 
(2S 3 + 1)(2S 4 + 1) p}' 



er 3+ 4->i + 2 (s) =0-1+2^3+4(8) 7^ , , i\ T' ( 14 ) 



and the role of the backward reactions ((cc)i+meson formation by D + D flavor 
exchange) can be explored without introducing any additional parameter once |Mo| 2 
is fixed. In Eq. (fTl"]) the quantities Sj denote the spins of the particles, while and 
p'j denote the cms momentum squared in the initial and final channels, respectively. 
The uncertainty in the cross sections (fT4"]) is of the same orde r of mag nitude as 
that in Lagrangian approaches using e.g. SU(4) flavor symm etry I ^ 44 | 14 5[ smcc t ne 
form factors at the vertices are essentially unknown EH . The cross sections for 
these backward channels - summed up again over all possible binary final states 
- are displayed in Fig. [20] separately for the 'non-strangeness' (upper part) and 
'strangeness' channels (lower part). 

The regeneration of charmonia by recombination of D (D*) meso ns in the 
hadronic phase was first studied by CM. Ko and collaborators in Ref. ^^ l The con- 
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Fig. 22. Schematic representation of a Pb + Pb collision at 160 A-GeV in space-time. The full 
dots represent early hard collision events (for Drell-Yan and cc-pairs) while mesons (it, n, p, lu, etc. 

- arrows) only appear after a respective formation time tp. The overlap area (inner rect angle) 
specifies the space-time region of hard production events. The figure is taken from Ref. 133 

elusion at that time was that this process was unlikely at RHIC energiesDS^EMIHI 
On the other hand, it was shown within HSD S3 that the contribution of the D + D 
annihilation to the produced J at RHIC is considerable. Moreover, the equilib- 
rium in the reaction J/ty + m <-> DD is reached (i.e. the charmonium recreation is 
comparable with the dissociation b y 'com oving' mesons). The reason for such dif- 
ferences is that the pioneering study E21 within the hadron gas model was confined 
to J/^f reactions with pions only and into two particular DD channels (D + D* 
and D* + D*). As one can see in Fig. [21j t he cr oss sections used in Refs. | 132 | 147 | 

- as obtained in the quark exchange model Q^Q] - in the two dissociation channels 
J/* + 7r — > D + D* and J/* + 7r — > D* + D* agree with the parametrization ([TO]). 
However, in HSD the interactions with all mesons into all possible combinations of 
DD states have been taken into account (cf. Fig. [T8|) . Note that the p-meson density 
at RHIC is large such that the cha nnel with the most abundant p-meson resonance 
is dominant. Furthermore, in Ref. ^\ the feed down from Xc and decays has 
been co nsidered. The res ults of Ref. S3 are in accordance with independent studies 
in Refs. H4 $ I 149 | 150 | 151 | s t rC ss the importance for DD annihilation in the late 
(purely hadronic) stages of the collisions. 

Note that in the default HSD approach (i.e. in the hadronic comovcr disso- 
ciation and recombination scenario) only formed comoving mesons participate in 
dissociation or DD recombination reactions (cf. Subsection l6.3p . Being hard probes, 
cc pairs are created in the early stage of the collision, while the comoving mesons 
are formed at a later stage. This situation is illustrated in Fig. [22] for a central 
Pb+Pb collision at 160 A-GeV for freely streaming baryons (thick black and red 
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Fig. 23. The calculated rate of J/ 'J dissociation reactions with mesons (solid histogram) for central 
Pb + Pb collisions at y/s = 17.3 GeV in comparison to the rate of backward reac tions of open 
charm pairs to J/4" + meson (dashed histogram). The figure is taken from Ref. 1^3 

lines). The initial string formation space-time points are indicated by the full dots; 
the mesons (indicated by arrows) hadronize after a time delay tp ps 0.8 fm/c as 
shown by the first hyperbola. A cc-pair produced in the initial hard nuclcon-nuclcon 
collision cannot be absorbed by mesons in the (green) shaded areas in space and 
time; however, a sizeable fraction of cc pairs (which should be produced within the 
inner rectangles) can also be produced in a dense mesonic environment. The upper 
hyperbolas in Fig. [2U represent the boundaries for the appearance of mesons from 
the second interaction points (full dots) which appear somewhat later in time; they 
stand for a representative further nucleon-nucleon collision during the reaction. 

Thus, the interaction with comoving mesons can occur only in the late stages 
of the reaction, i. e. after tp in their rest frame. This is taken into account in HSD 
by treating 'formed' particles and 'pre-hadrons' differently: 

• A 'formed' meson (baryon) is a quark-antiquark (quark-diquark) correla- 
tion - produced at time t$ - with hadronic quantum numbers that satisfies 
the two constrains: 

the time since production t — to > tp = jrp 
and 

the local energy density e < 1 GeV/fm 3 , 

where Tp 0.8 fm/c is the formation time, and 7 the Lorcntz 7-factor 
w.r.t. to the calculational frame. 
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Fig. 24. The calculated rate of J/* dissociation reactions with mesons (solid histogram) for central 
Au + Au collisions at y/s = 200 GcV in comparison to the rate of backward reac tions of open 
charm pairs to J/4< + meson (dashed histogram). The figure is taken from Ref. S3 

• A 'pre-hadron' is a correlation with hadronic quantum numbers for 

t — to < tF = 'yTF 

or 

e > lGcV/fm 3 . 

• leading quarks (diquarks) of the strings do not carry hadronic quantum 
numbers; they interact with reduced cross sections, which arc set to 1/3 
(2/3) of the hadronic ones in line with constituent quark number scaling. 

The result for the total J/vP comover absorption rate (solid histogram) in central 
Pb + Pb collisions at 160 A GeV is shown in Fig. [23] in comparison to the J/^S 
reformation rate (dashed histogram) that includes all backward channels. Since the 
rates differ by about 2 orders of magnitude, the backward rate for J/ 1 ^ formation 
can be neglected at SPS energies even for central Pb + Pb reactions. For central 
Au + Au collisions at = 200 GeV, however, the multiplicity of open charm pairs 
should be ~ 16, i.e. by about 2 orders of magnitude larger, such tha t a m uch higher 
J/* reformation rate (~ N^ s ) is expected at RHIC energies (cf. Ref.^SZJ). I n Fig. [24] 
we display the total J comover absorption rate (solid histogram) in comparison 
to the J/vI/ reformation rate (dashed histogram) as a function of time in the center- 
of-mass frame. Contrary to Fig. [23] now the two rates become comparable for t > 
4-5 fm/c and suggest that in central collisions at the full RHIC energy of ^/s ~ 
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200 GeV the J/vP comover dissociation is no longer important since the charmonia 
dissociated in this channel are approximately recreated in the backward channels. 

6.2. 'Threshold melting' (and recombination) 

The 'threshold melting' scenario is based on the idea of a sequential dissociation 
of charmonia with increasing temperature I25l 26 | 27 | 28 l j_ e Q f charmonium melting 
in the QGP due to color screening as soon as the fireball temperature reaches the 
dissociation temperatures (« 2T C for J/^f, ~ 1.1 — 1.2T C for excited states, where T c 
stands for the critical temperature of the deconfinement phase transition) . In the 
geomet rical Gl auber model of Blaizot et al. ^^ as well as in the percolation model 
of Satz I 27 | 154 |^ ^ j s assumCQ ; that the QGP suppression sets in rather abruptly as 
soon as the energy density exceeds a threshold value e c , which is a free parameter. 
These models are motivated by the idea that the charmonium dissociation rate is 
drastically larger in a quark-gluon-plasma (QGP) than in a hadronic medium E3 
such that further (hadronic) comover absorption channels might be neglected. 

We modify the standard sequential dissociation model of Refs. I 153 | 27 | 1 54] m 
two aspects: (i) the energy density is calculated locally and microscopically instead 
of using schematic estimates (c/. Section [3|); (ii) the model incorporates a char- 
monium regeneration mechanism (by DD annihilation processes as described in 
Section 6.1). The 'threshold scenario' for charmonium dissociation is implemented 
in HSD in a straight forward way: whenever the local energy density s(x) is above 
a threshold value Ej, where the index j stands for J/^, x c , the charmonium is 
fully dissociated to c + c. The default threshold energy densities adopted are 

ej/m = 16 GcV/fm 3 , e Xc = 2 GcV/fm 3 , and e 9 > = 2 GeV/fm 3 . (15) 

The dis sociation of charmon ia has been widely studied using lattice QCD (1QCD) 
in Refs. I 155 | 156 | 157 | 158 | 1 59| m orc [ er to determine the dissociation temperature (or 
energy density) via the maximum entropy method. On the other hand one may use 
potential models - reproducing the charmonium excitation spectrum in vacuum - 
to calculate Mott transition temperatures in a hot medium. Both approaches have 
their limitations and the quantitative agreement between the different groups is 
still unsatisfactory: 

• (A) Potential models usually employ the static heavy quark-antiquark 

pair free energy F(T) - calculated on the lattice - to obtain the char- 

moni um spect ral functions. This leads to the (low) dissociation tempera- 
tures [l6Qlfe] 

T me i t (J/«£) < 1.2 T c , T melt { Xc ) < T c , T meU (*') < T c . 

An alternative way is to use the internal energy U = F + TS as a 
quark-antiquark potential, which due to large contributions from the 
entropy S provides dissociation temperatures closer to the estimate 
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m gB2|163|164|165|16H17| 

It is presently unclear (and very much de- 
bated) how to extract a proper quark-antiquark potential from 1QCD cal- 
culations. 

• (B) The maximum entropy method is used to relate the Euclidean thermal 
correlators of charmonia - calculated on the lattice - to the cor respo nding 
spectral functions. This yields higher dissociation temperatures 

T melt (J/V) = 1.7-2 T c , T melt ( X c) = 1.1-1.2T C 

or [HSl 

T melt (J/V) > 1.5 T c , Uc)«i.ir c . 

Since the low values for the melting temperatures from the potential models are 
already in conflict with the J/^f data at SPS, the values IT5|) are employed in the 
following (if not stated otherwise). 



6.3. Charm interactions with pre-hadrons in the early phase 

In the default HSD approach all newly produced hadrons (by string fragmentation) 
have a formation time of Tp w 0.8 fm/c ps 1/Aqcd in their rest frame and do 
not interact during the 'partonic' propagation. Furthermore, hadronization is in- 
hibited, if the energy density - in the local rest frame - is above 1 GeV/fm 3 , which 
roughly corresponds to the energy density for QGP formation in equilibrium. This 
default approach underestimates the elliptic flow of hadrons at RHI C energies E^l, 
the suppression of hadrons wi th hi gh transverse momentum px ESI as we ll as the 
suppression of the far-side jets in central Au + Au collisions at y/s = 200 GeV. 
This failure has been attributed to the lack of explicit partonic interactions in the 
early collision phase, which corresponds to the phase of high energy density with 
the majority of hadrons being still 'under formation'. In order to simulate partonic 
interaction effects the HSD approach has been extended by ex plic it interactions of 
pre-hadrons with the (perturbative) charm degrees of freedom l^l 

Accordingly, an additional scenario has been implemented in the HSD simu- 
lations which is closely related to the 'comover suppression' scenario outlined in 
Scction [6-H which clearly separates 'formed hadrons' (existing only at energy den- 
sities below the energy density e cut = e c w 1 GeV/fm 3 ) from possible pre-hadronic 
states at higher energy densities, i.e, above the parton/hadron phase transition. 
Indeed, it is currently not clear whether D- or _D*-mesons survive at energy densi- 
ties above e c , but hadronic correlators with the quan tum numbers of the hadronic 
states are likely to persist above the phase transition I 170 [ Thus 'comovers' - with 
modified spectral functions - could show up also at energy densities above s c . We 
recall that in HSD a pre-hadron is defined as a state with the quantum numbers of a 
hadron, if the local energy density is above r* 1 GeV/fm 3 or if the state is still under 
'formation', i.e. if the time between production and hadronization is smaller than 
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the formation time rp (in its rest frame). For a more detailed d escription of the 
pre-hadron concept we refer the reader to S ection 16. II and to Rcfs. HZ 3- | 168 | 172 | 169 [ 
In line with the investigations in Refs. QSHESTthc J/ty production and ab- 
sorption in Au + Au collisions at y/s = 200 AGeV has been studied assuming 
the abso rpt ion of charmonia on pre-hadrons as well as their regeneration by pre- 



loni 

hadrons This adds additional interactions of the particles with charm quarks 
(antiquarks) in the very early phase of the nucleus- nucleus collisions as compared to 
the default HSD approach. Since these pre-hadronic (color-dipole) states represent 
some new degrees-of-frccdom, the interactions of charmed states with these objects 
have to be specified separately 

For notation, we define a pre-hadronic state consisting of a quark-antiquark pair 
as pre-meson rh and a state consisting of a diquark-quark pair as pre-baryon B. The 
dissociation cross section of a cc color dipole state with a pre-baryon is taken to be 
of the same order as with a formed baryon, 

= 5.8 mb, (16) 

w hereas th e cross section with a pre-meson follows from the additive quark model 

as mmm 

Elastic cross sections are taken as 



o% = 1.9 mb, o% m = (18) 



2 

Furthermore, elastic interactions of a charm quark (antiquark) are modeled by the 
scattering of an unformed D or D* meson on pre-hadrons with only light quarks as 



0^ = 3.9 mb, < n = \y e DB- 



In this way one can incorporate in HSD some dynamics of quark-antiquark pairs 
with a medium that has not yet formed the ordinary hadrons. However, it has to 
be stressed that further explicit partonic degrees of freedom, i.e. gluons and their 
mutual interactions as well as gluon interactions with quarks and antiquarks, have 
not been taken into account explicitly so far. 



6.4. Thermal and statistical models 

The statistical and thermal models are based on the assumption of statistical equi- 
librium, where the physical system can be characterized by a few Lagrange pa- 
rameters that specify the average energy as well as the average particle number, 
flavor content, etc. Since the fireball created in relativistic nucleus-nucleus collisions 
is rapidly expanding, equilibrium rate equations for charmonium dissoc iatio n and 
regeneration arc folded over a thermally evolving background in time In a 
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simplified form, the rate equation for the time evolution of a charmonium state ip, 
N^(r) can be written as 



dN, 



dr 



— — r^> 



Ni,~N^ . (20) 



The first key quantity in (|20|) is the inelastic charmon ium reaction rate, T^, which, 
by detailed balance, governs both gain and loss terms This is in quite analogy 
to the covariant transport models. In the QGP, the leading-order (LO) process is 
the well-kno wn gl uo-dissociation, g + — > c + c. However, as has been first empha- 
sized in Ref. H52l the gluo-dissociation process becomes inefficient for small J/ tp 
binding energies as expected due to color screening in the QGP (and even with- 
out screening for tp' and Xc states). Therefore, the quasi -free dissociation process, 
p + i\) — > c + c + p (p = q, q, g) has been introduced which naively is of next- 
to-leading order in a s , but provides a much larger phase space, and, consequently, 
the dominant dissociation rate for small charmonium binding (for gluo-dissociation, 
the phase space vanishes in the limit of vanishing binding energy). The other key 
quantity is the charmonium equilibrium limit, 7V^ 9 (r), which depends on the charm 
content and temperature of the system. The typical procedure is to assume charm 
production to be a hard process and thus to be restricted to primordial N-N colli- 
sions. The statistical model is then used to distribute the fixed number of cc pairs 
over the available charmed states in the system. This introduces both temperature 
and volume dependencies into N^ q (T(r)), as well as a sensitivity to medium mod- 
ifications of the charm states (e.g., red uced £)- meson masses lead to a reduction in 
the charmonium equilibrium numbers) | 122 | 174 | Alternative absorption mechanisms 
might al so play a role, su ch as gluon scattering on color dipole states as suggested 
in Refs. | 125 | 152 | 175 | 176 | or charmonium dissociation in the strong color fields of 
overlapping strings ^^3. 

The statistical hadronization model Q^Q] (originally proposed in Ref. E3) fol- 
lows a very different idea. The cc are produced by initially hard NN scattering, 
but charmonium bound states are generated i n a statistical fashion at the phase 
boundary from the QGP to the hadronic phase E3. This is fundamentally different 
from a regeneration of charmonia in the hadronic phase by D + D collisions and 
implies that all charmonia are emerging from a QGP phase. Since this assumption 
is questionable for peripheral nucleus-nucleus collisions, the authors have later on 
introduced 'corona effects' that mimic non-QGP charmonium formation. We will 
come back to actual results and comparisons in Section 8. 

Slightly later Th ews e t al. developed a coalescence model for charmonium pro- 
duction in the QGP H2SI w hich states that the abundance of charmonia increases 
roughly with the cc density squared and accordingly becomes increasingly impor- 
tant at higher bombarding energy (especially at LHC energies). While actual data 
at RHIC do not strongly support this scenario, it will be of substantial interest if 
the future dat a at LHC will provide evidence for charmonium enhancement instead 
of suppression | 1^0 | 173 [ 
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Fig. 25. The J/^f suppression as a function of the transverse energy Et in Pb+Pb collis ions at 
160 A-GeV. The solid line shows the HSD res ult w ithin the comover absorption scenario HZBB] 
The different symbols stand for the NA 50 da ta I I from the year 2000 (ana lysis A,B,C) while the 
dashed histogram is the UrQMD result The figure is taken from Ref. ESI. 

7. Anomalous suppression of J/ty: Comparison to data 
7.1. J/ty and suppression at SPS energies 

As pointed out before, 'cold nuclear matter' absorption effects alone cannot 
reproduce the strong suppression of J/\& observed by the NA50 collaboration in 
central Pb+Pb collisions (c/. Fig fTTj) . The extra suppression of charmonia in the high 
density phase of nucleus-nucleus collisions then may be attributed to the different 
'hot matter' scenarios outlined in Section [6l 

As a reminder, we recall the early J/^P suppression results for Pb+Pb at 160 
A- GeV (in the comover suppression scenario) both from the UrQMD and HSD 
transport calculations, which arc in line with the data of the NA50 Collaboration 
as demonstrated in Fig. [25j where the calculated J suppression is shown as 
a function of the transverse energy Et- The s olid line stands for the HSD result 
(within the comover absorption scenario) I 47 | 48 |^ wn ji e the various data points reflect 
the NA50 data from the year 2000 (anal ysis A,B,C) that agree reasonably well with 
the HSD and UrQMD calculations l^ 2 !! (dashed histo gram in Fig. 1251) . 

The anomalous suppression obs erved in S+U and Pb+Pb collisions by the 
NA38 ^ and NA50 Collaborations I 4 | 128 | j-^g Deen experimentally confirmed by 
NA60^in In+In collisions at 158 A-GeV. A couple of models have predicted J/^ 
suppression in In+In collisions as a function of centrality at 158 A-GeV based on 
the parameters fixed for Pb+Pb reactions at the same bombarding energy. How- 
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NA60,In+In, 158 A GeV 
QGP threshold melting 
NA50, Pb+Pb, 158 A GeV 
QGP threshold melting 




NA60, In+In, 158 A GeV 
Comover absorption 
NA50, Pb+Pb, 158 A GeV 
Comover absorption 




-100 



Fig. 26. The ratio S M(J cr( J/^t) /cr(DY) as a function of the number of participants N par t in In+In 
(red line with open squares) and Pb+Pb reactions (blue line with open circles) at 158 A- GeV 
relative to the normal nuclear absorption given by the straight black line. The full dots and squares 
denote the respective data from the NA50 and NA60 Collaborations. The model calculations reflect 
the comover absorption model (right part) and the 'QGP threshold scenario' (left part) with 
e J/* = 16 GcV/fm 3 , e Xc = 2 GeV/fm 3 , e^/ = 2 GeV/fm 3 while discarding comover absorption. 
The figure is taken from Ref . EH 



ever, the predictions within the 'Glaube r based' comover model and the 'threshold 
melting scenario' from Refs.^ ^^^^^^ have failed to describe the In+In data with 
sufficient accuracy (cf. Section 7.2). 

The charmonium production and suppression in In+In and Pb+Pb reactions 
at SPS energies has been reinvestigated in the HSD transport approach in 2006 
within the 'hadronic comover model' as well as the 'QGP threshold scenario' EH. 
As found in Rcf. EH the comover absorption model - with a single parameter |Mo| 2 
for the matrix element squared for charmonium-meson dissociation - performs best 
with respect to all data sets for J /iff suppression as well as for the iff' to J /if! 
ratio for Pb+Pb (cf. Figs. I26I27|) . We recall that the iff' suppression is presented 
experimentally by the ratio 

B^(J/V -h. w)a(J/*)/a(DYy 1 ' 

In the HSD calculations this ratio is taken as 0.0165 for nucleon-nucleon collisions, 
which is again based on the average over pp, pd,pA reactions^. The centrality 
dependence of the ratio is then a prediction of the model for different systems and 
bombarding energies. 

The 'QGP threshold scenario' roughly reproduces the J /if! suppression for both 
systems at 158 A-GeV (Fig.[26j l.h.s.) apart from the very central Pb+Pb collisions. 
The comover absorption model follows slightly better the fall of the J /if! survival 
probability with increasing centrality (Fig. [26j r.h.s.), whereas the 'threshold sce- 
nario' leads to an approximate plateau in both reactions for high centrality 

On the other hand, the 'threshold melting scenario' clearly fails for the iff' to 
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Fig. 27. The to J/\P ratio as a function of the transverse energy Et for Pb+Pb at 1 60 A-GeV. 
The full dots and stars denote the respective data from the NA50 Collaboration EE] The HSD 
resultHH for the comover absorption model is shown as the red line, whereas the blue line indicates 
the 'QGP threshold scenario' with Ej/y 
blue line with £j/m = 16 GeV/fm 3 , e Xc 
comover absorption). The figure is taken from Ref. 



16 GeV/fm 3 , e Xc = 2 GeV/fm 3 = Ey, and the light 
= 2 GeV/fm 3 and Ey, = 6.55 GeV/fm 3 (discarding 

m 



J/\I/ ratio, since too many \J/' already melt away for a critical energy density of 
2 GeV/fm 3 at 158 A-GcV (cf. Fig. [27j). Only when assuming the to dissolve 
above ~ 6.5 GeV/fm 3 a reasonable description of all data is achieved in the 'QGP 
threshold scenario'; this threshold, however, is not in accordance with lattice QCD 
calculations such that the 'threshold scenario' meets severe problems. 

Indeed, the extra suppression of charmonia by comovers is seen in Fig. [28] (solid 
red lines) to match the J/^! suppression in In+In and Pb+Pb as well as the to 
J/* ratio (for Pb+Pb) rather well. The more recent data (1998-2000) for the 
to J/\I r ratio agree with the HSD prediction SHI within error bars. This had been 
a problem in the past when comparing to the 1997 data (dark green stars). One 
may conclude that the comover absorption model so far cannot be ruled out on the 
basis of the available data sets from the SPS within error bars. The to J/^ ratio 
for In+In versus centrality is not yet available from the experimental side but the 
theoretical predictions (provided in Fig. ]28|) might be verified/falsified in future. 

Some com ment s on the comover absorption model appear in place: As shown in 
Fig. 7.2 of Ref. ^Sl, the comover densities in central Pb+Pb collisions at 158 A-GcV 
become quite large and almost reach 2 fm -3 in the maximum, which appears high 
for 'free' mesons with an eigenvolume of about 1 fm 3 . However, as mentioned before, 
the quasi-particle mesons considered here dynamically should not be identified with 
'free' meson states that show a long polarization tail in the vacuum. As known from 



August 11, 2008 12:47 WSPC/INSTRUCTION FILE Linnyk 



42 Linnyk, Bratkovskaya, Cassing 



40 



3 30 



Q 



20 



10 



50 



0.015 



0.010 



3 



0.005 



0.000 



HSD 



i — ,— ln+ln, l^o A LreV |- 




Pb+Pb, 158 A GeV 




^O - Comover absorption 
+ NA50 2004 

►Ti 1 






Yjpj^ Comover absorption 
HNLl ■ NA60 2005 




HSD 







100 

N 



150 



200 



100 



part 



200 

N 



300 



400 



part 



• Comover absorption 




HSD 



. * NA50 1997 
• NA50 1998-2000 
Comover absorption 



50 100 150 200 25 50 75 100 125 150 



Fig. 28. The ratio B ll/J ,cr(J/ty)/o-(DY) as a function of the number of participants in In+In (l.h.s.) 
and Pb+Pb reactions (r.h.s.) at 158 A-GeV. The full symbols denote the data from the NA50 
and NA60 Collaborations (from Refs. I 4 l 5 l 10 i3). The solid (red) lines show the HSD results for the 
comover absorption model with a matrix element squared |Mq| 2 = 0.18 fm 2 /GeV 2 . The lower 
parts of the figure show the HSD results in the same limit for the \P' to J/St ratio as a function 
of Npart (for In+In) or the transverse energy Et (for Pb+Pb). The vertical lines on the graphs 
reflect th e th eoretical uncertainty due to limited statistics of the calculations. The figure is taken 
from Ref.HH 



lattice QCD, the correlators for pions and p-mesons survive well above the critical 
temperature T c , such that 'dressed' mesons, i.e spectral densities with the quantum 
numbers of the pseudo-scalar and vector (isovector) modes, also show up at high 
energy density (similar to the J/\& discussed above H 29 | 130 | 18 0fy Such 'dressed' 
mesons are expected to have a more compact size in space. 
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Fig. 29. Same as Fig.[28]but for the 'QGP threshold scenario' with e J/9 = 16 GeV/fm 3 , e Xc = 2 
GeV/fm 3 = £\$ji while discarding comover absorption, i.e. for |Afo| 2 — 0. The figure is taken from 



Rcf. 



44 



The results for the 'threshold scenario' for J /if? as well as if?' are displayed in 
Fig. [2U in comparison to the same data for the thresholds = 16 GeV/fm 3 , 

e Xc = 2 GeV/fm 3 = while discarding any dissociation with comovers, i.e. |Mo| 2 
=0. In this scenario, the J /if? suppression is well described for In+In, but the sup- 
pression is slightly too weak for very central Pb+Pb reactions. This result emerges, 
since practically all \ c and if?' dissolve for N part > 100 in both systems whereas 
the J /if? itself survives at the energy densities reached in the collision. Since the 
nuclcon dissociation is a flat function of N part for central reactions, the total ab- 
sorption strength is flat, too. The deviations seen in Fig.[29]might indicate a partial 
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Fig. 30. NA50 datal 128 l 18 3on J/ip (left panel) and ip' (right panel) production in Pb(158 AGeV)- 
Pb collision s at the SPS, compared to solutions of a kinetic rate equation in a thermal fireball 
background s tart ing from initial yields subject to primordial nuclear absorption. The figure 
is taken from Ref. 1221. 



melting of the J/^f for N part > 250, which is not in line with current lattice QCD 
calculations claiming at least ej/^ > 5 GeV/fm 3 . In fact, a lower threshold of 5 
GeV/fm 3 (instead of 16 GeV/fm 3 ) for the J /if? has practically no effect on the 
results shown in Fig. [29] Furthermore, a threshold energy density of 2 GeV/fm 3 
for the \I>' leads to a dramatic reduction of the 'J' to J/^f? ratio, which is in severe 
conflict with the data (lower part of Fig. [29]) . Also note that due to energy density 
fluctuations in reactions with fixed N part (or Ex) there is no step in the sup press ion 
of J /if? versus centrality as pointed out before by Gorenstein et al. in Ref. H&U. 



7.2. Thermal and statistical models at SPS 

It is of interest to compare the results of dynamical (transport) models with the 
statistical and thermal fireball models, which assume statistical equilibrium during 
the nucleu s- nuc leus collision. 

In Ref. 022 the rate equation (20]) has been solved for P6(158 AGeV)-P6 col- 
lisions using a fireball evolution fitted to transport calculations. The authors also 
include primordial nuclear absorption and suppression in the hadronic phase. The 
resulting centrality dependence for J/-0 production (including feed-down form Xc 
and ty') is shown in the left panel of Fig. 1301 The only free parameter in this ap- 
proach is the strong coupling constant entering into the quasi-free dissociati on cr oss 
section in the QGP which has been fixed to a s ~ 0.25. The NA50 data^Hl are 
fairly well reproduced, with a small contribution from regeneration (incorporated 
also by detailed balance). The main effect for the direct J/i/Ss is their suppres- 
sion which is largely restricted to the QGP (after nuclear absorption as inferred 
from p-A data). On the other hand, suppressi on is found to be substantially 
affected by the hadronic phase, and the NA50 data can on ]y b c reproduced in 
this approach, if in-medium (dropping) D-meson masses are implemented, which 
accelerate the direct if?' — ► DD decays, cf. right panel of Fig. [30] 
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Fig. 31. NA60 data for J/ip production in In(158 AGeV)-In ESI compared to theoretical model 
predictions that are in app roxim ate agreement with the NA50 Pb — Pb data: perc olation model 
(upper (dashed) blue line) 154 ^ kinetic rate equation (middle magenta line) EEH and hadronic 
comovers (lower green line) 11221 Data and theory curves are normalized to an "expected yield" 
which includes the effects of primordial nuclear absorption as extracted from p-A data. The figure 
is taken from Ref. UM\. 



A widely debated issue, furthermore, is whether the NA50 data support the no- 
tion of a more or less sharp "onset behavior" of J/ip suppression, possibly related 
to the formation of a deconfined medium. The NA60 collaboration has scrutinized 
this issue by measuring the J/i p sup pression pattern in a medium-size system, 
i.e., In-(158 AGeV)-In. The data I I, normalized to the J/ip yield expected after 
primordial nuclear absorption, are compared to theoretica l pre dictions in Fig. 1311 
dissociation of cc states in QGP by Satz, Digal, Fortunato ^^ (b lue l ines); regen- 
eration of charmonium in the QGP by Rapp, Grandchamp, 

BrownQS! 

(pink lines); 

comover abs orpti on by mesons in a Glauber type model (green lines) by Capella 
and Ferreiro The latter approaches have been adjusted to the NA50 Pb+Pb 
data, but it turns out that non e of them fully describes the In-In measurements. 
The hadronic comover scen ario 1 over-predicts the suppression throughout, the 
schematic percolati on m odel misses the onset significantly, while the kinetic rate 
equation approach somewhat over-predicts the suppression for the most cen- 
tral collisions. Overall, the predictions of the kinetic approach (in a time-dependent 
thermal fireball background) do not fare too badly with the data. 
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Fig. 32. The J/\E> nuclear modification factor Raa for Au+Au collisions at yfs = 20 G eV 
as a function of the number of participants N par t in comparison to the data from Ref. 12ZI for 
midrapidity (full circles) and forward rapidity (full triangles). HSD results for the 'QGP threshold 
melting' scenarios are displayed in terms of the lower (green solid) lines for midrapidity J/^'s 
(\y\ < 0.35) and in terms of the upper (orange dashed) lines for forward rapidity (1.2 < y < 
2.2) without recombination (l.h.s.) and with a recombination via D + D. The error bars on the 
theoretical results indicate the statistical uncertainty due to the finite number of events in the HSD 
calculations. Predictions for the ratio B^^^^a^i / -B Mfl (J j as a function of the number of 
participants N par t for Au+Au at ^fs = 200 GeV are shown in the lower set of plots. The figure 
is taken from Ref. EH. 



7.3. J/* and at RHIC 

Up to date a simultaneous description of the seemingly energy-independent sup- 
pression of J together with its narrow rapidity distribution and a strong elliptic 
flow v 2 of charmed hadrons - as found at RHIC - has presented a challenge to mi- 
croscopic theories. The large discrepancies of present studies are striking in view of 
the success of the hadron-string transport theories in describing charmonium data 
at SPS energies. This has lead to the conjecture that the sizeable difference between 
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the measured yields and transport predictions is due to a neglect of the transition 
from hadronic to partonic matter, i.e. the strongly-coupled Quark-Gluon-Plasma 



In the RHIC experiments, one defines the nuclear modification factor Raa as 



where dN A ' A jdy denotes the final yield of J/vP in AA collisions, dN pp jdy is the 
yield in elementary pp reactions while N co u is the number of initial binary collisions. 

In the upper part of Fig.[32]we present a comparison of Raa(J/^) (from HSD) 
for Au + Au c ollis ions as a function of the number of participants N par t to the 
data from Ref. ^3. The results for the 'threshold melting' scenario (without the 
reformation channels D + D — > ( J/^, Xc, + meson) are displayed on the l.h.s. 
of Fig.[32]in terms of the lower (green) solid line for mid-rapidity J/^'s (\y\ < 0.35) 
and in terms of the upper (orange) dashed line at forward rapidity (1.2 < \y\ < 2.2). 
The experimental data from PHENIX^^are given by the full circles at mid-rapidity 
and by triangles at forward rapidity. In this simple scenario practically all charmonia 
are dissolved for N part > 50 due to the high energy densities reached in the overlap 
zone of the collision, which is clearly not compatible with the PHENIX data and 
indicates that charmonium reformation channels are important. 

The reformation and dissociation channels (D + D <-> (J/fy,x c ,^') + meson) 
are switched on after a formation time. The results for this model study are dis- 
played in the upper right part of Fig. [3H and demonstrate that for N part > 200 
an approximate equilibrium between the reformation and dissociation channels is 
achieved. However, here the calculations for forward rapidity match the data at 
mid-rapidity and vice versa showing that the rapidity dependence is fully wrong. 
Furthermore, the J suppression at more peripheral reactions is severely over- 
estimated. Thus one has to conclude that the 'threshold melting + reformation 
scenario' is incompatible with the PHENIX data and has to be ruled out at top 
RHIC energies. 

In the lower parts of Fig. [32] the results for the ratio of the and J/^> dilcpton 
yields (given by their cross sections multiplied by the corresponding branching 
ratios) are shown, which so far have not been measured. Here the two recombination 
models give finite ratios as a function of centrality but predict a larger ty' to J /"P 
ratio at forward rapidity than at mid-rapidity which is a consequence of the higher 
comovcr density at mid-rapidity. Experimental data on this ratio should provide 
further independent information. 

The suppression of charmonia by the 'comover' dissociation channels within 
the model described in Section 16.11 is presented in Fig. I33( where the charmonium 
reformation channels by D + D annihilation have been incorporated. The HSD 
results for Raa in the purely hadronic 'comover' scenario are displayed in the upper 
part of Fig. [33] in comparison to the data from Ref. using the same assignment 
of the lines as in Fig. [32] This scenario gives a continuous decrease of Raa{ J /^) 
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Fig. 33. Same as Fig. 1321 for the 'comover abs orpt ion scenario' including the charmonium refor- 
mation channels. The figure is taken from Ref. 



with centrality. However, the rapidity dependence of the comover result is opposite 
to the one experimentally observed, as dictated by the higher comover density at 
mid-rapidity. 

The VP' to J/^f ratio is displayed in the lower parts of Fig. [33] and shows a 
decreasing ratio with centrality similar to the results at SPS energies (cf. Fig. |28|). 
As pointed in Ref. SH, an independent measurement of \&' will provide further 
information on the charm reaction dynamics and final charmonium formation. For 
instance, a leveling off of the to J ratio with increas ing centr ality would be a 
signal for charm chemical equilibration in the medium H3ZHZZIISS1 Additionally, it 
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Fig. 34. The ratio of the nuclear modification factors Raa at mid-rapidity (|y| < 0.35) and at 
forward rapidity (1.2 < \y\ < 2.2) vs centrality in Au + Au collisions at y/s = 200 GeV. The 
HSD results in the purely hadronic scenario ('comover absorption') are displayed in terms of the 
blue dashed line (with open circles) and in case of the 'threshold melting' scenario in terms of 
the violet dot-dashed line (with open squares). The error bars on the theoretical results indicate 
the statistical uncertainty due to the finite number of Monte-Carlo events in the calculations. The 
lower full green dots represent the data of the PHENIX Collaboration I2ZI Note that the data 
have an additional systematic uncertainty of ±14%. The lower solid (red) line with stars gives the 
result for the 'comover absorption' scenario when i nclu ding additional prc-hadronic interactions 
with charm (see text). The figure is taken from Ref . 1221 



provides a very clear distinction between the 'threshold melting' scenario and the 
'comover' approach. 

The non-applicability of the traditional 'comover absorption' model and 'thresh- 
old melting' scenario at the top RHIC energy is most clearly seen in the centrality 
dependence of the ratio of the nuclear modification factors Raa at forward rapid- 
ity (1.2 < \y\ < 2.2) and at mid-rapidity (|y| < 0.35) as shown in Fig. [3l The 
HSD results in the purely hadronic scenario ('comover absorption') are displayed 
in terms of the blue dashed line (with open circles) and in case of the 'threshold 
melting' scenario in terms of the dot-dashed magenta line (with open squares). The 
error bars on the theoretical results indicate the statistical uncertainty due to the 
finite number of Monte-Carlo events in the calculations. The lower full green dots 
in Fig. l34l represent the corresponding data of the PHENIX Collaboration which 
show a fully different pattern as a function of centrality (here given in terms of the 
number of participants N part ). The failure of these 'standard' suppression models 
at RHIC has lead to the conclusion in Ref. EH that the hadronic 'comover absorp- 
tion and recombination' model is falsified by the PHENIX data and that strong 
interactions in the pre-hadronic (or partonic) phase is necessary in order to explain 
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Au+Au, s =200 GeV 
Prehadron interactions 




Fig. 35. The J/9 nuclear modification factor Raa J22l for Au+Au collisions at = 200 AGeV 
as a function of the number of participants N par t in comparison to the data from I2ZI for mid- 
rapidity (full circles) and forward rapidity (full triangles). The HSD results for the hadronic 
'comover' scenario including additionally pre-hadronic interactions of charm according to 11161 - 
H91 arc displayed in terms of the upper (green solid) line with open circles for mid-rapidity J/9 1 s 
(\y\ < 0.35) and in terms of the lower (orange d ashe d^ line with open triangles for forward rapidity 
(1-2 < \y\ < 2.2). The figure is taken from Ref.ESl. 



the large suppression at forw ard rapidities. 

Consequently, in Ref. additional interactions of charm with pre-hadrons (as 
described in Section 6.3) have been incorporated in the 'comover scenario' to have 
a first glance at the dominant effects. The J/^ suppression pattern in this case 
is shown in Fig. [35] in comparison to the same data as in Fig. [33] for the hadronic 
comover model. When including the pre-hadronic interactions, the suppression pat- 
tern for central and forward rapidities becomes rather similar to the data within the 
statistical accuracy of the calculations. Indeed, the ratio of Raa at forward rapidity 
to mid-rapidity now follows closely the experimental trend as seen in Fig. [33] by the 
lower red solid line. 

Some further information may be gained from the J /\& rapidity distributions in 
Au+Au collisions at RHIC. The latter distribution is shown in Fig. [36] in compar- 
ison to the PHENIX data for central collisions (upper l.h.s.), semi-central (upper 
r.h.s.), semi-peripheral (lower l.h.s.) and peripheral reactions (lower r.h.s.) for the 
standard 'comover' scenario (dashed blue lines) and the 'comover' model including 
additionally pre-hadronic interactions of charm according to (fl"6j) - (fT9|) (solid red 
lines). Whereas for peripheral reactions these additional early interactions practi- 
cally play no role, the additional prc-hadron elastic scattering lead to a narrowing 
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Fig. 36. The rapidity distribution dNj/y /dy for different centralities from the standard 'comover' 
model (dashed blue lines) and the 'comover' model with additional prc-hadronic interactions of 
charm according to 11161 - 11191 (solid red lines). The full dots show the respective data from the 
PHENIX Collaboration E3. The calculated lines have been smoothed by a spline algorithm. The 
reactions are Au+Au at y/s = 200 GeV. The figure is taken from Ref. l^2J 



of the J rapidity distribution with the centrality of the collision (roughly in line 
with the data). In the standard 'comover' model an opposite trend is seen: here the 
interactions of charmonia with formed hadrons produce a dip in the rapidity distri- 
bution at y « which increases with centrality since the density of formed hadrons 
increases accordingly around mid-rapidity. Since the total number of produced cc 
pairs is the same (for the respective centrality class) and detailed balance is incor- 
porated in the reaction rates, we find a surplus of J at more forward rapidities. 
The net result is a broadening of the J/^> rapidity distribution with centrality in 
the purely hadronic scenario opposite to the trend observed in experiment. 

In summarizing the results at RHIC energies, the hadronic 'comover' dynamics 
for charmonium dissociation and recreation as well as the charmonium 'melting' 
scenario do not match the general dependence of the J/^> on rapidity and central- 
ity as seen by the PHENIX Collaboration. In fact, a narrowing of the J/\& rapid- 
ity distribution cannot be achieved by comover interactions with formed hadrons, 
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since the latter appear too late in the collision dynamics. Only when including 
early pre-hadronic interactions with charm, a dynamical narrowing of the charmo- 
nium rapidity distribution with ccntrality can be achieved, as demonstrated within 
the pre-hadronic interaction model. Consequently, the PHENIX data on J /if? sup- 
pression indicate the presence and important impact of pre-hadronic or partonic 
interact ions in the early charm dynamics. This finding is in line with earlier studies 
in Refs. |47| 168 | 169 | demonstrating the necessity of non-hadronic degrees of freedom 
in the early reaction phase of nucleus-nucleus collisions at RHIC energies for the 
description of the elliptic flow v 2 , the suppression of hadrons at high transverse 
momentum pr and far-side jet suppression. 

7.4. Charmonia excitation functions 

In this Section we present the excitation functions for the J/if? survival probability 
in Au + Au collisions from FAIR to top RHIC energies in the different scenarios in 
order to allow for a further distinction between the different concepts. The results of 
HSD calculations are presented in the upper part of Fig. [37] for the 'QGP threshold 
melting + hadronic recombination' scenario (dashed green lines with open triangles) 
and the 'comover absorption + recombination' model (solid red lines with open 
circles) for central (l.h.s.) and minimum bias (r.h.s.) Au+Au reactions as a function 
of the beam energy. The J/if? survival probability Sj/^ is defined as 



where N^ n and N B ' B denote the final number of J /if? mesons and the number 
of J/if?'s produced initially by BB reactions, respectively. We find that from FAIR 
energies of 20 - 40 A- GcV up to top SPS energies of 158 A-GcV there is no 
significant difference between the two models for the J /if? survival probability in 
case of central collisions. The differences here show up mainly in the full RHIC 
energy range where the 'QGP threshold melting + hadronic recombination' scenario 
leads to substantially lower J/if? survival probabilities. In case of minimum bias 
collisions the 'comover absorption + recombination' model (solid lines) gives a 
roughly energy independent J/if? survival probability, whereas the 'QGP threshold 
melting + hadronic recombination' scenario shows lower J/if? survival probabilities 
(lower dashed green lines) for laboratory energies above ~ 100 A-GeV due to a 
larger initial melting of J/vP-mesons at high energy density. 

A clearer distinction between the different concepts is offered by the excitation 
functions for the if?' to J/if? ratio in Au + Au collisions. The calculated results 
are shown in the lower part of Fig. [37] for the 'QGP threshold melting + hadronic 
recombination' scenario (dashed green lines with open triangles) and the 'comover 
absorption + recombination' model (solid red lines) for central (l.h.s.) and mini- 
mum bias reactions (r.h.s.). Here the iff' is already melting away in central Au+Au 
reactions in the 'QGP threshold melting' scenario at bombarding energies above 40 
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Fig. 37. upper part: The excitation function for the J/ty survival probability in the 'QGP threshold 
melting + hadronic recombination' scenario (dashed green lines with triangles) and the 'comover 
absorption + recombination' model (solid red lines with circles) for central (l.h.s.) and minimum 
bias Au+Au reactions (r.h.s.) as a function of the beam energy. Lower part: The VP' to J/VB 
ratio for the same reactions as in the upper part of the figure in the 'QGP threshold melting + 
hadronic recombination' scenario (dashed green lines with triangles) and the 'comover abs orpt ion 
+ recombination' model (upper solid red lines with circles ). The figure is taken from Ref. 



A-GeV, whereas a substantial amount of 'J' survives in the 'comover absorption + 
recombination' model. Thus measurements of VP' suppression at the lower SPS or 
top FAIR energies will clearly distinguish between the different model concepts. 
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Fig. 38. The survival probability Sj/ip (left plot) and ratio \&' to J/^S (right plot) as a function 
of the number of participants N par t in Au+Au reactions at 25 A-GeV. The blue lines (with 
open dots) reflect the 'threshold scenario' for £j/<j = 16 GeV/fm 3 , e Xc = 2 GeV/fm 3 , e^r = 
6.55 GeV/fm 3 while the violet line (the lower line with open dots on the r.h.s.) stands for the 
'threshold scenario' for = 16 GeV/fm 3 , e Xc = 2 GeV/fm 3 , e^i = 2 GeV/fm 3 . The solid red 

lines (full dots) denote the results for the comover absorption model with the standard matrix 
element squared |Mo| 2 = 0.18 fm 2 /GeV 2 . The dashed line (l.h.s.) represents the HSD calculations 
including only dissociation channels with nucleons. The figure is taken from Ref. HH. 



7.5. FAIR energies 

The CBM Collaboration at GSI is aim ing a t charmonium measurements in heavy- 
ion collisions at the future FAIR facility ^^ . This opens up the possibility to explore 
the charmonium suppression mechanism at lower bombarding energies of about 25 
A- GeV in Au+Au collisions. The corresponding HSD predictions are displayed in 
Fig. [38] for the survival probability Sj/q, (l.h.s.) and the ratio ^' to J/^ (r.h.s.) as 
a function of the number of participants N part . The violet line in Fig. 1381 stands for 
the 'threshold scenario' with ey, i.e. e,//^ = 16 GeV/fm 3 , e Xc = 2 GeV/fm 3 , e^,/ 
= 2 GeV/fm 3 . The solid red lines denote the results for the comover absorption 
model with the standard matrix element squared |Afo| 2 = 0.18 fm 2 /GcV 2 . 

We note that in Au+Au reactions at 25 A-GeV in the 'threshold scenario' (solid 
line with open dots in the left plot) only a very low amount of Xc and no J 
are melted at the energy densities reached in these reactions. On the other hand 
the comover density decreases only moderately when stepping down in energy from 
158 A- GeV to 25 A- GeV such that the J/\I/ survival probability in the comover 
absorption model (lower solid line in the left part) is lower. This is even more 
pronounced for the 'J' to J /^S ratio versus centrality, which in the 'threshold melting 
scenario' (middle line in the right part) is clearly above the result achieved in 
the comover absorption model (lower line in the right part). Consequently, the 
different dissociation scenarios may well be distinguished in future charmonium 
measurements at FAIR. 
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Fig. 39. Ratio of the averaged J/* to n multiplicity for Pb + Pb at the SPS beam energy of 158 
A-GeV at mid-rapidity (l.h.s.) and in full 4tt acceptance (r.h.s.) as a function of the number of 
binary collisions N co u for the different suppression scenarios implemented in HSD - the 'comover' 
model (dashed blue line with open circles) and the 'threshold melting' scenario (green dot-da shed 
line with open triangles) - in comparison to the statistical model by Gorenstein and GazdzickiES] 
(r.h.s.; straight orange line) and the statistical had ronization model by Andronic et al. HSS1 (l.h.s.; 
solid black line) . The figure is taken from Ref . 1221 



8. Testing the assumption of statistical hadronization 

The assumption of statistical hadronization - i.e. of J/^'s being dominantly pro- 
duced at hadronization in a purely statistically fashion according to available phase 
space and the number of available c and c quarks - leads to a scaling of the 
(J/^)/(h) ratio with the system size 033 where (h) is the average hadron mul- 
tiplicity. Since (h) ~ (ir), the ratio (J/^)/(n) has been calculated in HSD in the 
different scenarios for charmonium suppression: 

• 'threshold melting' + recombination via DD — > cc + m including the back- 
ward reactions cc + m — > DD, 

• hadronic ('comover') absorption: DD — > cc+m and the backward reactions 
cc + m — > DD; 

• 'prchadron interactions': DD — > cc + m and the backward reactions cc + 
m — > DD as well as early prc-hadronic charm interactions as described in 
Section 6.3. 

The results of the calculations are shown in Fig. [39] toge ther with the prediction 
of the statistical model of Gorenstein and Gazdzicki for the full phase space 
(st raight or ange line; r.h.s.) and the statistical hadronization model by Andronic et 
al. I 186 | 188 | for m id-rapidity (solid black line; l.h.s.) for Pb+Pb at 158 A-GcV. The 
centrality dependence here is given by the number of initial binary collisions N co u. 
The actua l comparison in Fig. [39] indicates that the statistical model by Andronic 
et al. E?6] predicts a sizeably larger J/ty to it ratio at midrapidity for peripheral 
and semi-peripheral reactions than the microscopic HSD results for the different 
scenarios. For central reactions - where an approximate equilibrium is achieved - 
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Fig. 40. Same as Fig. [39] but for Au + Au at the top RHIC energy of ^ = 200 GeV. The red 
solid line shows additionally the result of the 'comov er' model including the pre-hadronic charm 
interactions (see text). The figure is taken from Ref. 



all scenarios give roughly the same ratio. In full A-k phase space the HSD results 
indicate also a slightly higher J/ty to it ratio in the 'comover' model relative to the 
'melting' scenario but both ratios only weakly depend on centrality - roughly in line 
with the statistical model of Gorenstein and Gazdzicki IEH (orange straight line). 
Consequently, only peripheral reactions of heavy nuclei might be used to disentangle 
the different scenarios at top SPS energies at midrapidity (or in full phase space). 

The situation is different for Au+Au collisions at the top RHIC energy as may 
be extracted from Fig. [4171 where the J/^S to pion ratio (l.h.s.: at mid-rapidity; r.h.s.: 
for 47T acceptance) is shown as a function of N co u. The standard 'comover' model 
(dashed blue lines) is only presented for reference but is unrealistic according to 
the analysis in Section 7. We find that the 'comover' model with early pre-hadronic 
charm interactions ( solid red line with stars, l.h.s.) is very close to the statistical 
hadronization model (solid black line) at mid-rapidity except for very peripheral 
collisions. The 'threshold melting' scenario follows the trend in centrality but is 
down by about 30%. Thus at mid-rapidity there is no essential extra potential in 
differentiating the scenarios. Considering the full 4tt acceptance (r.h.s.) one finds a 
practically constant J/\& to pion ratio for N co u > 200 from the HSD calculations 
as expected fr om t he statistical model, however, the early model of Gorenstein 
and Gazdzicki is down by about a factor of ~10 (and may be ruled out by 
present dat a) . The latter conclusion is in agreement with the independent analysis 
in Ref. 033. 



9. Transverse mass spectra 

Apart from rapidity dependent particle abundances, also the dynamics in the trans- 
verse direction (to the beam) provides relevant information. 
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Fig. 41. Upper plot: The transverse mass spectra from pp collisions at Ti a f, = 25 GeV for pions (full 
squares), kaons (open triangles), and 0-mesons (full rhombi) from the string model as implemented 
in HSD. The dashed line shows an exponential with slope parameter Eq = 0.143 GeV. Middle 
plot: The same as the upper plot, but for pp reactions at T; a j, = 160 GeV. The dashed line shows 
an exponential with slope parameter Eq = 0.176 GeV. Lower plot: The same as the upper plot, 
but for pp reactions at y/s = 200 GeV. The dashed lin e shows an exponential with slope parameter 
E = 0.225 GeV. The figures are taken from Ref.EH 
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9.1. Elementary collisions 

We recall again that HSD successfully reproduces the measured transverse mass 
= \/p\ + m? x spectra of pions and kaons in pp collisions for bombarding en- 
ergies of 12 GcV, 24 GeV as well as for collisions at ^ = 200 GeV^. In Fig. gT] 
we present additionally the corresponding spectra for </>-mesons, D + D mesons and 
charmonia at SPS, RHIC and FAIR energies. We display the differential multiplici- 
ties (2rriT)~ 1 dNx / drriT in the transverse mass. The pion spectra describe the sum 
of 7r + , 7T°, 7r~, the kaon spectra the sum of K + , K°, K°, K~ , the Z?-meson spectra 
the sum of all D, D* , D s , D* and their antiparticles while the spectrum denoted by 
cc includes the J/^, the Xc as well as the where the latter contribution starts 
at niT ~ 3.7 GeV and becomes visible as a tiny kink in the my-spectra. Here the 
open charm and charmonia results stem from the parametrizations specified in Sec- 
tion [2] (including the decay \c — ► J/*& + 7), while the spectra for pions, kaons and 
^-mesons are from the LUND string model as implemented in the HSD transport 
approach. For orientation, we also show exponential spectra with slope parameters 
of 143 MeV, 176 MeV and 225 MeV, respectively, which describe the my-spectra 
of pions rather well. The kaon spectra at all energies are down by a factor of ~ 3, 
the (j> spectra by a factor of 9-10 relative to this line due to strangeness suppres- 
sion in pp collisions. However, it is quite remarkable that the charmonia spectra 
fit well to this approximate TOy-scaling (within a factor of 2-3) at yfs = 7.1, 17.3 
and 200 GeV, respectively. Furthermore, the spectrum of open charm is roughly 
compatible with my-scaling at y/s = 17.3 and 200 GeV, while the D,D mesons 
are suppressed relative to the scaling by a factor ~ 30 close to threshold (y/s = 7.1 
GeV). Such an 'apparent' statistical pro duct ion of mesons in elementary reactions 
has been advocated before by Becattini HS2I. We have to stress, however, that all 
these observations on the charm sector are based on our extrapolations (Section ^ 
and have to be checked experimentally. 

9.2. SPS and RHIC energies 

As mentioned before (cf. Ref.^3) the hadron spectra at high are underestimated 
both in HSD and UrQMD compared to the data at energies from SPS to RHIC. 
Nevertheless, it is illustrative to have a global view on the transverse mass spectra 
of hadrons with different flavor in case of central Au+Au collisions. In this respect 
Fig. [42] displays the results of the hadronic transport at SPS and RHIC energies 
which illustrates an approximate meson 7«T-scaling for all flavors (cf . Rcf. . The 
mr-scaling for pions, kaons, D-mesons and J/\V in central collisions of Au + Au at 
top SPS and RHIC energies is essentially due to an approximate my-scaling in pp 
collisions at yfs = 17.3 GeV and substantial D, D and J/\V final state interactions 
in the nucleus-nucleus case. 

We recall that the measured transverse mass spectra of hadrons (heavier than 
pions) at AGS, SPS and RHIC energies show a 'hardening' in central Au+Au colli- 
sions relative to pp interactions (cf. Ref. for low transverse mass or momentum. 
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Fig. 42. Left panel: The transverse mass spectra of pions (full squares), kaons (open triangles), 
0-mesons (full rhombi), D + D mesons (open squares) and J/ty,^' mesons (full dots) in the HSD 
approach for a central Au + Au collision at 160 A-GeV. The thin dashed line shows an exponential 
with slope parameter Eq = 0.176 GeV. Right panel: Same as left panel, but at top RHIC energy 
of 21.3 A-TeV. The thin dashed l ine shows an exponential with slope parameter Eq = 0.225 GeV. 
Both figures are taken from Ref. \32\. 



This hardening is commonly attributed to collective flow, which is absent in the 
respective pp or pA collisions. Consequently it is important to get precise data on 
open charm and charmonium transverse momentum (j>t) spectra, since their slope 
migh t give information on the pressure generated in a possible early partonic phase 
H2QI. This argument is expected to hold especially for J/^> mesons, since their elastic 
rescat tering cr oss section with hadrons should be small in the hadronic expansion 
phase I 191 I 19 S1 ( see a i so Section HOI). 



9.3. FAIR energies 

In contrast to the observation at SPS energy, the approximate mr-scaling for pions, 
kaons, Z?-mesons and J/vI/ no longer holds for central collisions of Au+Au at 25 
A- GeV as demonstrated in Fig. [43l Here the HSD calculations show a suppression 
of D-mesons by a factor of ~ 10 relative to the global m-p-scaling - characterized 
by a slope of 143 McV - if no D-meson self energies are accounted for. 

On the other hand, attractive mass shifts of D, D mesons of -50 McV at po 
might be expected due to had ronic interaction models when extending SU(3) flavor 
to SU (4) flavor symmetry 11221. At the densities of 5-8 po chiral symmetry should be 
restored, i.e. the large < qq > condensate of the nonperturbativc vacuum should 
have disappeared. Accordingly, the production of cc pairs in the 'new' perturbativc 
vacuum might be enhanced since only the invariant mass of a cc has to be produced 
e.g. by gluon-gluon fusion. Accordingly, the threshold for DD production - which 
is ~ 3.739 GeV in vacuum - might be reduced by w 2 ■ 0.35 GeV = 0.7 GeV in the 
chirally restored phase to about 3 GeV, only. Such a reduction of the cc production 
threshold leads to an enhancement of open charm mesons by about a factor of 7 
(magenta crosses in Fig. I43j) such that an approximate m^-scaling for all mesons 



August 11, 2008 12:47 WSPC/INSTRUCTION FILE Linnyk 



60 Linnyk, Bratkovskaya, Cassing 




m T [GeV] 

Fig. 43. The transverse mass spectra of pions (full squares), kaons (open triangles), <j>- mesons (full 
rhombi),D + D mesons (open squares) and J/\&, W mesons (full dots) in the HSD approach for a 
central Au+Au collision at 25 A- GeV without including self energies for the mesons. The crosses 
stand for the D-meson mj spectra when including an attractive mass shift of — 50p/po MeV. The 
thin dashed line shows an exponential with slope parameter Eo = 0.143 GeV. Note that final state 
elastic scattering of kao ns a nd i^-mesons with pions has been discarded in the calculations. The 
figure is taken from Ref. 1221 

is regained. Thus, a global scaling of all mesons may be regarded as a strong 
medium effect on the charmed hadrons and as a signature for a chirally restored 
phase - if observed by CBM. 

10. High pt quenching of open charm and charmonia 

A significant suppression of high transverse momentum hadrons in Au+Au colli- 
sions compared to pp is observed at RHIC energies of y/s = 200 GeVl 1 ^^^^^ 
and is attributed to the ene rgy loss of highly energetic particles in a hot col- 
ored mediu m (QGP)EQll 2Q2l i n f act; the recent observation by the PHENIX^31 
STAR ^231 and BRAHMS collaborations that a similar suppression is not ob- 
served in d+Au interactions at mid-rapidities at the same energy supports this 
idea. 

Vigorous theoretical efforts are under way to understand parton energy loss in 
terms of perturbativc QCD (pQCD). Various groups have described the suppression 
of light hadrons in terms of radiative energy loss by gluon bremsstrahlung. Accord- 
ing to such calculations, charm and beauty quarks should be absorbed significantly 
less than light quarks and gluons due to their higher mass. However, data from the 
PHENIX and STAR experiments, which compare the production in nucleus-nucleus 
to proton-proton collisions of high-p^ "non-photonic" electrons (which originate 



August 11, 2008 12:47 WSPC/INSTRUCTION FILE Linnyk 



Charm in heavy ion collisions 61 



1/2 

D+Dbar from Au+Au, s =200 GeV, mid-rapidity 



4 PHENIX ('07) 0-10% centrality 
■ PHENIX ('06) 0-10% 
□ 10-20% 
* STAR ('05) 0-5% 



syst +- 



PHENIX ('06), 20-40% 
•HSD,b=5fm 




syst+/-8% 



PHENIX ('06), 40-60% 
— C — HSD, b=7 fm 



PHENIX ('06), 60-92% 
■HSD,b=12fm 



0.5 . 





syst +/- 12% 



syst+/-28% 



1 2 3 4 5 
P T [GeV/c] 



1 2 3 4 5 
P T [GeV/c] 



Fig. 44. HSD predictions for the ratio of the final to the initial (i.e. at the production point) 
transverse momentum spectra of D + S-mesons (solid lines with open dots, color: blue) from 
Au + A u collisions at y/s = 200 Ge V for b = 1, 5, 7 and 12 fm at mi d-rap idity as calculated in 
Ref.lSH The PHENIX data from Ref. HSU (denoted '06) and from Ref. IMil (denoted '07) on R AA 
of non-photonic electrons as well as the STAR data from Ref. ESS have been added later. 



mainly from heavy-flavor decays) seem to indicate that heavy quarks loose energy 
in a comparable fashion as light quarks. 

In order to quantify the effect of hadronic final state interactions, we show in 
Fig. 23] the HSD predictions from Ref. f or the ratio of the final to the initial 
transverse pt spectra of D +D- mesons (solid lines with open dots, color: blue) from 
Au + Au collisions at ^fs = 200 GeV calculated for impact parameter b — 1,5,7 and 
12 fm at midrapidity. HSD predicts an enhancement of D,D mesons at low momenta 
with a maximum at pr ~ 1 GeV/c and a relative suppression for px > 2 GeV/c. 
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Fig. 45. HSD predictions for the ratio of the final to the initial pt spectra of J/^-mcsons (solid 
lines with open dots, color: red) from A u + Au collisions at y/s = 200 GeV for b = 1, 5, 7 and 12 
fm a t m id-rapidity as calculated in Ref. The preliminary PHENIX data for Raa(J/^) from 
Ref. I2ZI have been added later. The grey dashed bands represent the results from the extended 
version of the HSD comover approach, in which prchadron interactions are taken into account as 
described in Section 16.31 



These effects increase with the centrality of the Au+Au collision. We note that the 
maxima in the ratios disappear when switching off the rescattering with mesons in 
the transport approach. Thus a collective acceleration of the D + D mesons occurs 
also via elastic scattering with mesons. As seen in Fig. |44]the suppression seen by 
PHENIX may well be explained by hadronic comover interactions up to transverse 
momenta about 4 GeV/c. Only for higher pt a clear signal for parton energy loss 
- cither gluon bremsstrahlung or parton elastic scattering - may be extracted in 
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comparison to the data! 

The suppression pattern Raa( J /^) from HSD (in the comover scenario) is quite 
analogous to that of Z?-mcsons for different ccntrality showing a slight maximum for 
transverse momenta of ~ 2 GeV/c and a steady decrease for higher pr- In Fig. [45] 
the predictions (from Ref. for the ratio of the final to the initial transverse px 
spectra of J/^-mesons (solid red lines with open dots) as well as new calculations 
in the extended version of the comover approach - as described in Section [6.3l fgrev 
dashed bands) - from Au + Au collisions at s/s = 200 GcV for four different 
centrality regions at mid-rapidity are displayed. The preliminary PHENIX data 
for Raa{J/^) from Ref. - added later - show a substantially different pattern 
especially for non-peripheral interactions. The strong suppression for low px J/^! 
mesons seen experimentally suggests that not primordial J/'J's are accelerated 
during the dynamical evolution but that at least a part of initially formed J/^'s 
are dissolved and created later e.g. by cc coalescence. We stress that the reformation 
of charmonia in the hadronic phase (by D + D etc.) carries the flow from the D- 
mesons and thus does not lead to suppression at small px as seen experimentally. 
This observation supports the idea that part of the charmonia are produced in the 
hadronization process! 



11. Collective flow 

A further possible way to disentangle hadronic from partonic dynamics is the elliptic 
flow 1*2(2/, Pt) as a function of the rapidity y and transverse momentum px- The flow 
^2(2/1 Pt) is driven by different pressure gradients in case of non- vanishing spatial 
anisotropy 

e 2 =< 4^4 > • (24) 

Since €2 decreases fast during the expansion of a nonccntral nucleus-nucleus re- 
action, the magnitude of gives information about the interaction strength or 
interaction rate of the early medium. 

The phenomenon of collective flow can generally be characterized in terms of 
anisotropics of the azimuthal emission pattern, expressed in terms of a Fourier series 

dN 

— (if) oc 1 + 2v! cos(0) + 2v 2 cos(20) + . . . (25) 
dip 

which allows a transparent interpretation of the coefficients Vi and v-x- The dipolc 
term v\ arises from a collective sideward deflection of the particles in the reaction 
plane and characterizes the transverse flow in the reaction plane (the "bounce-off" ) . 
The second harmonic describes the emission pattern perpendicular to the reaction 
plane. For negative V2 one has a preferential out-of-plane emissio n, called squeeze- 
out. Pions at SIS energies exhibit a clear out-of-plane preference I 205 | 206 | w ] 1 i c i 1 [ s 
due to shadowing by spectator nucleons^l 
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Fig. 46. Elliptic flow V2 of J/'J's produced in central and peripheral In+In collisions at 158 A-GeV 
beam energy in the hadronic 'comover' mode of HSD (open circle and open triang le - ) c ompared to 



the NA60 data 1211] represented by black diamonds. The figure is taken from Ref. 
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P resen tly, the most employed flow observables are the in-plane and elliptic 
flows ^ 



vi = < — 



Pi 

PT 



v — ( — ^ y \ (26) 
2 \Pl +P 2 v J ' 

Here, p x denotes the momentum in a;-direction, i.e. the transverse momentum within 
the reaction plane and p y the transverse momentum out of the reaction plane. The 
total transverse momentum is given as px 



lp x + Py] the z-axis is in the beam 
direction. The bounce-off, the squeeze-out and the 

antiflow [20l20i( third 

flow com- 



ponent I2**f) have been suggested as differential barometers for the properties of com- 
pressed, dense matter from SIS to RHIC. In particular, it has been shown E-0 | 208 | 
that the disappearance or "collapse" of flow might be a direct result of a first order 
phase transition. 



11.1. SPS energies 

In Fig. |3S]we compare the HSD result for ^(J/^) at SPS in the purely hadronic 
'comover' scenario in co mparison to the data for V2(J/^) of the NA60 collaboration 
for In+In collisions EH j n central collisions the elliptic flow is practically zero both 
in the calculation as well as in the experiment, whereas in peripheral reactions a 
nonzero flow emerges. The agreement (within error bars) between the theory and 
the data indicates that, in line with the reproduction of the J/* suppression data^ 
(see Section l7.1j) . the low amount of V2 does not point towards additional strong 
partonic interactions at SPS energies. Consequently, the present measurements of 
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Fig. 47. Elliptic flow of D-mesons produced in Au + Au collisions at = 200 GeV as a f unct ion 
of pt from HSD (solid blue line with open circles) in comparison to the PHENIX data CSS on 
V2 of non-photonic electrons. The red line with open stars shows the HSD result for the V2 of D- 
mesons when including additionally prc-hadronic charm interactions as described in Section 16.31 
The figure is taken from Ref . 1221 

J/\E' elliptic flow (at SPS energies) do not provide further constraints on the model 
assumptions. 

11.2. RHIC energies 

The situation, however, is different for the collective flow of D-mesons at top RHIC 
energies. Though collectiv e flow can be described very elegantly in hydrodynamics 
(c/. Refs.l21i ? | 216 | 217 | 218 |) by a proper choice of initial conditions, one has to be very 
careful, since most hydrodynamical calculations - describing flow - fail to reproduce 
the hadron spectra with the same initial condition s (an d vice versa). Hydrodynamic 
flow and shock formation has been proposed early^^ as the key mechanism for the 
creation of hot and dense matter during relativistic heavy-ion collisions. However, 
the full three-dimensional hydrodynamic al flow problem is much more complicated 
than the one-dimensional Landau model 12121 u §ed in many of the present hydro- 
dynamical calculations. The 3-dimensional compression and expansion dynamics 
yields complex triple differential cross-sections, which provide quite accurate spec- 
troscopic handles on the equation of state. In this respect, it is important to consider 
also microscopic multi-component (pre-) hadron transport theory (see Section[2]) as 
control models for viscous hydro and as background models to extract interesting 
non-hadronic effects from data. 

In Fig.|47]we show the elliptic flow of D-mesons produced in Au + Au collisions 
at y 7 ! = 200 GeV as a function of the transverse momentum pr in HSD (solid blue 
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Fig. 48. Elliptic flow of J/'I'-mesons produced in Au + Au collisions at y/s = 200 GeV as a function 
of Pt from HSD. The red line with open stars shows the HSD result for the V2 of Z?-mesons when 
including additionally pre-hadronic charm interactions as described in Section 16.31 
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Fig. 49. Elliptic flow of J/\l/-mesons produced in Au + Au collisions at \fs = 200 GeV as a function 
of pt from HSD and other m odels . Preliminary PHENIX dataEEEEEHl are shown by symbols. 
The figure is taken from Ref . ESI. 



line with open circles) compared to the PHENIX data on V2 of non-photonic 
electrons. Here the elliptic flow of D-mesons is clearly underestimated in the dc- 
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fault (purely hadronic) HSD model (c/. Ref. l^). Only when including pre-hadronic 
charm interactions - as described in Section 6.3 - the elliptic flow moderately in- 
creases (red line with open stars) , but still stays below the PHENIX data. We thus 
conclude that the modeling of charm interactions by pre-hadronic interactions - 
as described in Section 6.3 - accounts for part of the non-hadronic generation of 
the V2, but does not provide enough interaction strength in the early phase of the 
collision. Quite remarkably, this find ing i s again fully in line with the underestima- 
tion of high pt hadron suppression USE] a s well as far-side jet suppression USS] m 
the pre-hadronic interaction model. Independently, also the charm collective flow 
points towards strong partonic interactions in the early reaction phase beyond the 
pre-hadronic scattering incorporated so far. 

Since a large fraction of J/^'s in central Au+Au collisions at RHIC are cre- 
ated by D — D recombination, the elliptic flow of J/vP's obtained from HSD in 
the comover (purely hadro nic) case is c omparatively small, too (cf. Fig. |48|) . The 
preliminary PHENIX data | 212 | 213 |214| additionally shown by blue symbols in 
Fig. 221 where the HSD prediction is represented by a brown (dash-dot) line and a 
band. The accuracy of the preliminary data so far does not allow for a differentiation 
between the different model predictions presented in Fig. H§l 

In principle, partonic cascade simulations l ^0 | 221 | 222 | 223 | 224 | 225 | 226 | ^ 
promising tools in approaching the mechanism of the early generation of elliptic 
flow. However, unexpectedly high parton cross sections of ~ 5-10 mb have to be 
assumed in parton cascades ED m order to reproduce the elliptic flow seen ex- 
perimentally. These cross sections (per constituent quark) are roughly the same 
in the partonic and hadronic phase (or even higher in the parto nic phase). In the 
new version of the parton cascade model of Xu and Greiner 1222] lower binary cross 
sections are employed but additional 2 3 gluonic reactions are incorporated that 
are very sensitive to the screening masses (parameters) used. Furthermore, a strong 
elliptic flow of particles may als o be generated by repulsive mean-fields as known 
from the nuclear physics context ^2l. In all the parton cascades addressed above no 
mean- fields or parton selfcnergics arc incorporated. In fact, a recent study within 
the Parton-Hadron-String-Dynamics (PHSD) approach suggests that a large frac- 
tion of the final elliptic flow might stem from strong repulsive partonic mean-fields 



11.3. FAIR energies 

The expected results for the _D-mcson elliptic flow V2 at mid-rapidity for Au+Au 
reactions at 25 A-GeV from HSD is compared in Fig. [50] with the i>2 of charged 
hadrons. As seen from 1501 the D, D elliptic flow is smaller than the V2 of the lighter 
hadrons following essentially the mass ordering. Such a low elliptic flow is due to 
the moderate interaction rates of D, D mesons in the hadronic expansion phase. A 
measurement of a sizeably larger elliptic flow of Z?-mesons at FAIR would indicate 
the presence of partonic degrees of freedom already at bombarding energies of 25 - 



43 
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Fig. 50. The HSD predictions for the elliptic flow v% of D + S-mesons (solid lines with full dots) 
and charged hadrons (solid lines with open triangles) from Au + Au collisions at 25 A-GeV at 
b = 7 fm versus pt for \y\ < 1. The figure is taken from Ref . 12QI 

35 A-GcV. 

12. Conclusions 

This review has been devoted to an update of the information on charm dynamics 
gained from p+A and A+A reactions at SPS and RHIC energies in comparison to 
microscopic transport studies by incorporating different scenarios for the anomalous 
suppression of charmonia seen experimentally in central heavy-ion reactions. We 
have found that present data on charmonium suppression for Pb+Pb and In+In 
reactions at top SPS energies (158 A-GeV) compare well with transport calculations 
in the hadronic comovcr model involving only a single parameter for the average 
matrix element squared |Af | 2 that fixes the strength of the charmonium cross 
sections with comoving hadrons. This holds for the J /"P suppression as well as the 
iff' to J/\& ratio versus collision centrality. The bare 'QGP threshold scenario' gives 
satisfying results for the J /if? suppression for both systems (In+In and Pb+Pb) at 
158 A-GcV but fails in the *f?' to J/^f? ratio since too many if?' already melt away for 
a critical energy density of 2 GeV/fm 3 at 158 A-GeV. These findings suggest that 
the charmonium dynamics in heavy-ion reactions is dominantly driven by hadronic 
interactions in the SPS energy regime. Since energy densities above 1 GeV/fm 3 are 
reached in central nucleus-nucleus collisions at 158 A-GeV, our observation indicates 
that hadronic correlators (with quantum numbers of the familiar hadrons) still 
persist above the critical energy density for the formation of a QGP. This finding 
is independently supported by lattice QCD calculations that show a survival of 
strongly bound hadronic states up to temperatures of 2 T c , which corresponds to 
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an energy density of about 30 GeV/fm 3 . In fact, such high energies densities (or 
temperatures) are not reached even in central Au+Au (or Pb+Pb) collisions at the 
SPS. 

The situation is found to be essentially different at top RHIC energies of y/s = 
200 GeV (or 21.3 A-TeV). The study of the formation and suppression dynamics of 
J r /\&, Xc and ^' mesons within the HSD transport approach for Au+Au reactions 
at the top RHIC energy has demonstrated that both scenarios i.e. 'charmonium 
melting' in a QGP state as well as the hadronic 'comover absorption and recre- 
ation model' fail severely at RHIC energies. This is found in the J/^S suppression 
pattern versus the centrality of the collision, in the J/^f rapidity distribution as 
well as in the differential elliptic flow of J/^B and the charmonium nuclear mod- 
ification factor Raa as a function of transve rse momentum px (in comparison to 
recent data from the PHENIX Collaboration^^). Especially the latter observable 
indicates that at least part of the final J/\I/'s are created by coalescence of cc pairs 
in the hadronization phase. Only when including prc-hadronic degrees in the early 
charm reaction dynamics, the general suppression pattern of charmonia may be 
reasonably described; though, the elliptic flow V2 is still (slightly) underestimated. 
On the other hand, Raa(j>t) for <//\& mesons cannot be described in the comover 
approach, even when incorporating the early prchadron interactions. This analysis 
demonstrates that the dynamics of c, c quarks in heavy-ion reactions at RHIC en- 
ergies are dominated by partonic interactions in the strong QGP (sQGP) and can 
neither be modeled by 'hadronic' interactions nor described appropriately by color 
screening alone. 

The different scenarios of 'charmonium melting' and 'hadronic comover disso- 
ciation' may clearly be distinguished at FAIR energies (of about 25 - 35 A-GcV), 
where the centrality dependence of the J /iff survival probability and the to J /"P 
ratio are significantly lower in the 'comover absorption' model than in the 'charmo- 
nium melting' scenario. This result comes about since the average comover density 
decreases only moderately with lower bombarding energy, whereas the region in 
space-time with energy densities above critical values of e.g. 2 GeV/fm 3 decreases 
rapidly and ceases to exist below about 20 A-GcV even in central collisions of 
Au+Au. 

13. Open questions and perspectives 

Although substantial progress has been achieved on the field of charmonium physics 
in the last years - both theoretically and experimentally - the description of charmo- 
nium and open charm dynamics in relativistic nucleus-nucleus collisions remains a 
challenging task. So far we have just obtained a very rough picture of the collective 
flow and attenuation pattern of charmed and hidden-charm mesons and a profound 
microscopic understanding is still lacking. 

On the theoretical side the basic open problem - and future challenge - is to 
incorporate explicit partonic degrees of freedom in the description of relativistic 
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nucleus-nucleus collisions and their transition to hadronic states in a microscopic 
transport approach. The available string/hadron transport models only schemati- 
cally incorporate some parton dynamics (by modeling interactions of prc-hadrons) 
but definitely fail in describing the partonic phase in accordance with an equation 
of state from lattice QCD. On the other side the present partonic cascade sim- 
ulations (propagating massless partons) fail in describing the reaction dynamics 
when employing cross sections from perturbative QCD. One might argue that this 
is due to the neglect of partonic sclfcncrgics which even in case of small interaction 
strength should be large due to the high parton densities reached in ultra-relativistic 
nucleus- nucleu s collisio ns. In fact, the studies of Pcshicr in a dynamical quasipar- 
ticle approach I 228 | 229 | mc ji ca t e that the effective degrees of freedom in a partonic 
phase should be quite massive and have a width in the order of the pole mass al- 
ready slightly above T c . Such short lived degrees of freedom have to be propagated 
in off-shell transport approaches whi ch also allow for a description of hadronization 
in terms of 'local' rate equations S3. 

On the experimental side, further differential spectra of charmonia and open 
charm mesons will become available from RHIC-II and the LHC in the next years. 
While studies at RHIC most likely probe the 'partonic fluid' produced in the heavy- 
ion reactions, experiments at LHC - especially for high pt degrees of freedom - might 
provide further surprising results with respect to the properties of the partonic 
system at extreme energy densities. Also the future charm measurements at FAIR 
are expected to provide new insight into the charm dynamics especially at high 
baryon density. Accordingly, the field of open charm and charmonium physics is far 
from being closed and will remain exciting for a long time. 
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